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ABSTRACT 



The United States Coast Guard area and district commands are 
assigned platforms to accomplish required missions. As the nature 
and profile of the missions change over time, the platforms must be 
reallocated to better satisfy the new mission needs. The problem 
of reallocating existing platforms and, perhaps, acquiring new ones 
is referred to as the force structure problem. This problem is 
complex because of the multi-mission capability of the platforms, 
and the difficulty in quantifying projected mission requirements 
and platform suitability in meeting these requirements. 

In this thesis, methods for quantifying the mission 
requirements and platform suitability are proposed. Using these 
methods, a linear integer optimization model to allocate platforms 
is developed. Additionally, a FORTRAN based interface is 
implemented to facilitate the utilization of the model. Included 
in this interface are basic functions of a database system to aid 
the users in maintaining and updating model data. To demonstrate 
the use of the system, a realistic sample problem was extracted 
from Fiscal Year 1989 operational and administrative data. Outputs 
for the problem are given and analyzed. 
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THESIS DISCLAIMER 

The reader is cautioned that computer programs developed in this 
research may not have been exercised for all cases of interest. 
While every effort has been made, within the time available, to 
ensure that the programs are free of computational and logical 
errors, they cannot be considered validated. Any application of 
these programs without additional verification is at the risk of 
the user. 
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I. 



INTRODUCTION 



The origin of the United States Coast Guard's motto SEMPER 
PARATUS . Always Ready , was lost in antiquity. The Coast Guard 
now faces many political, technological, and economic 
challenges that must be addressed to ensure that this motto 
remains an accurate description of the Coast Guard. The 
proper management of capital assets is necessary for the well- 
being of any organization. The Coast Guard's Force Structure 
problem can be broadly described as striving for the efficient 
procurement and operation of the Coast Guard's current and 
future fleet of cutters 1 and aircraft (or, more generally, 
platforms) to accomplish the myriad of mandated programs and 
missions. 

A. BACKGROUND 

The United States Coast Guard (USCG) is both operationally 
and administratively decentralized with its headquarters 
located in Washington DC. Two area commands, in New York and 
San Francisco, direct operations in the Atlantic and Pacific 
regions respectively. Area Commanders maintain operational 
control over cutters greater than 180 feet in length [Ref. 1] . 



1 Cutters is the Coast Guard term for describing vessels 
sixty-five feet or longer. 
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The area commands are further divided into ten separate 
geographical commands called districts (see Figure 1) . All 
aircraft and smaller cutters are allocated or assigned to 
districts. Each district is controlled by a Flag Officer who 
exercises autonomous authority over the platforms assigned to 
the district and uses them to accomplish missions within the 
district's boundaries. 




Coast Guard duties and responsibilities are historically 

organized into seven major operating programs Defense 

Readiness, Search and Rescue (SAR) , Enforcement of Laws and 
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Treaties (ELT) , Marine Environmental Protection (MEP) , Aids to 
Navigation (ATON) , Ice Operations (IOP) , and Marine Safety. 
Each major operating program consists of one or more missions. 
Coast Guard missions are defined as "major subdivisions of 
Coast Guard statutory [operating programs] which are detailed 
in U.S. Law" [Ref. 2:p. 2]. To obtain "Coast Guard resource 
activity data for use in facility planning/management, program 
cost allocation and program management," the missions are 
further subdivided and expanded into approximately twenty 
employment categories [Ref. 3]. Employment categories are 
used to define the specific uses of a platform, including both 
training and operations. Table 1 provides an example of the 
relationship between major operating programs, missions, and 
employment categories as defined by the Coast Guard. 

Although confusing and technically incorrect, it is common 
practice to refer to major operating programs as missions. 
Perhaps, this is due to the fact that the term mission is more 
readily identified with activities to be accomplished and, in 
essence, major operating programs are general groupings of 
these activities. This thesis will deal with requirements 
computed at the major operating program level. To avoid 
confusion, this thesis will differentiate programs and 
missions as defined in the previous paragraph or in Table 1. 
The abbreviations assigned in Table 1 will be used during this 
thesis and in the sample problem discussed in Chapter V. 
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TABLE 1. 

MAJOR OPERATING PROGRAMS, MISSIONS AND EMPLOYMENT CATEGORIES 



Major Operating 
Programs 


Missions 


Employment 

Categories 


Marine 

Environmental 

Protection 

(MEP) 


Marine 

Environmental 

Response 


Marine 

Environmental 
Response | 

Marine Inspection 


Defense Readiness 
(MILOP) 


Port Safety and 
Security 


Port Safety and 
Security 


Defense Readiness 


Military Ops 

Military Training 

Military 

Exercises 


Enforcement of 
Laws and Treaties 
(ELT) 


Fisheries 

Enforcement 


Domestic 

Fisheries 

Foreign Fisheries 


Drug Interdiction 


Drugs - Air 
Interdiction 

Drugs - Surface 
Interdiction 


Enforcement of 
Laws and Treaties 


Immigration Laws 
ELT - Other 


Aids to 

Navigation 

(ATON) 


Short Range Aids 
to Navigation 


Short Range Aids 
to Navigation 


Radio Navigation 


Radio Navigation 


Ice Operations 
(IOP) 


Polar Operations 


Polar Operations 


Domestic Ice 


Domestic Ice 


Marine Safety 
(MSAFE) 


Commercial Vessel 
Safety 


Marine Inspection 
Program 


Boating Safety 


Boating Safety 


D Search and Rescue 
I (SAR) 


Search and Rescue 


Search and Rescue 
1 



4 



To accomplish all missions in the seven major operating 
programs, the Coast Guard owns and operates over 2,500 
platforms consisting of approximately 300 cutters and over 200 
fixed wing aircraft and helicopters. The remaining platforms 
are smaller vessels assigned to shore facilities (i.e., groups 
and stations) for performing near coastal and inland 
operations [Ref. 4:p. 150]. Most Coast Guard cutters and 
aircraft have a multi-mission capability to respond to "a wide 
variety of maritime activities and shift emphasis on short 
notice when need arises" [Ref. 5:p. 17]. This allows for 
greater flexibility in operations and procurement. However, 
this multi-mission capability of platforms also increases the 
complexity of assigning platforms to efficiently and 
effectively accomplish the missions in the seven operating 
programs on an annual basis. This problem of assigning 
platforms to missions, or directly to districts, is referred 
to as the force structure problem by the Coast Guard. 

B. RELATED AND PRIOR WORK 

In the past, the Coast Guard has attempted to quantify and 
solve the Force Structure problem, but lacked the necessary 
training or tools to analyze alternative solutions effectively 
[Ref. 6:p. 8]. The Coast Guard is currently without an 
acceptable analysis of its platform requirements (or force 
structure) , though several projects are either being developed 
or proposed to assist in solving this formidable problem. 
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Most of the ongoing projects currently focus on specific 
missions and/or platforms. The following sections provide a 
brief description of relevant studies and ongoing projects. 

1. Previous Coast Guard Studies 

Prior to 1989 the Coast Guard annually developed two 
separate planning documents, cutter and aviation requirements, 
addressing the force structure problem. The purpose of these 
documents was to be: 

...an analytical planning document which develops the 
Coast Guard ' s ... capital resources needed to achieve 
goals established in approved Operating Program Plans 
(OPP' s) . . . . [and] support formulation of the Coast Guard 
Capital Investment Projection, which assigns priorities 
over a five year time span [Ref. 7 & 8]. 

Though these two documents were independently 
prepared, they used similar methodologies. These 
methodologies are based on having individual program managers, 
one for each of the seven major operating programs, specify 
annual requirements in terms of employment hours for each 
district. To illustrate the concept of employment hours, 
consider a cutter which just completed one hour of an 
Enforcement of Laws and Treaties (ELT) mission. If the cutter 
monitors Search and Rescue (SAR) frequencies during part of 
the mission, the commanding officer may report that the cutter 
completed 1.5 employment hours, one hour for ELT and 0.5 hours 
for SAR, even though the actual mission time is only one hour. 

There are two sources of data used in determining the 
number of required platforms for each district. One is from 
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the seven program managers who specify the number of 
employment hours to be performed by each platform type in 
order to accomplish missions in their programs. The other set 
of data is the hours required for platforms in each district 
to perform support programs such as public affairs, training 
and research and development. These hours are extracted from 
the data in the previous three years of the Abstract of 
Operations 2 . Summing the hours from the two sources gives the 
total employment hours required for each platform type in each 
district. In an attempt to compute the actual usage hours, 
the employment hours are multiplied by predetermined factors 
known as Simultaneous Mission Factors (or SMF) (see, e.g. , Ref 
7 , page 2-4) . 

To determine the number of platforms required by each 
district, the Coast Guard simply divides the actual usage 
hours as computed above by the number of hours one platform 
can operate. If the number of required platforms of a 
particular type exceeds the number of platforms currently 
assigned to the district, then substitute platforms would be 
assigned to accomplish the hours instead. 

In October 1989, the Cutter and Aircraft Rebaseline 
Requirements Study (or the Rebaseline Study) was completed and 
intended to replace the cutter and aviation requirement 
documents. The Rebaseline Study used a method similar to the 

2 A quarterly operations statistics report submitted by 
all Coast Guard platforms [Ref. 3] 
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ones in these documents. The main difference is that the 
total employment hours are based on estimates prepared by the 
members of the Rebaseline Study group, district and area 
commanders as well as the program managers. 

Major drawbacks with the above studies are as follows. 
First, platforms selected or specified by program managers or 
the Rebaseline Study to perform a particular mission does not 
necessarily imply the efficient usage of the available 
platforms. This is due to the fact that some missions may be 
effectively performed by unselected platforms which may 
otherwise be under utilized. Second, specifying the number of 
mission employment hours required by each platform type in 
each district or area predetermines the number of platforms to 
be allocated to that district or area. Third, the concept of 
employment hours is subjective for it depends solely on the 
judgement of officers commanding the platforms [Ref. 3:p. 4- 
3]. Lastly, simultaneous mission factors which accurately 
convert employment hours to actual hours of usage are 
difficult to determine. These drawbacks may have undermined 
the acceptance of the Rebaseline Study. 

2. Ongoing Projects 

The Coast Guard has initiated two similar simulation 
projects to examine platform utilization in the performance of 
two major missions; SAR and ELT. Since June 1986, a KSS 
(knowledge-based decision support systems) project has been 
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under development to assist decision makers in the acquisition 
of cutters [Ref. 6]. Recently, a DSS (Decision Support 
System) was developed to examine the patrol boat (WPB) 3 Force 
Structure and Homeporting issues using simulation for the 
underlying models [Ref. 9]. The Coast Guard Research and 
Development (R&D) Center has requested Fiscal Years 92-93 
funding for a "General Force-Mix" project "founded upon an 
optimization model" [Ref. 10] . Through these projects, the 
Coast Guard has gained a better understanding of the Force 
Structure problem and is beginning to use the decision tools 
made available through mathematical programming and management 
science techniques. 

C. STUDY GOAL 

The primary goal of this thesis is to develop a 
methodology to assist the Coast Guard in determining an 
effective platform allocation on an annual basis. This 
methodology involves modeling the allocation process as an 
optimization model and quantifying the mission requirements 
for each district. 

It is hoped that the results described in this thesis 
illustrate the potential of mathematical programming models as 
decision aids and further encourage the Coast Guard to employ 
the methodology in the future. 



3 Patrol boats are Cutters 80 to 120 feet in length 
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D. THESIS ORGANIZATION 



Chapter II describes a method of quantifying annual 
mission requirements and platform utilization. Chapter III 
outlines the formulation of an integer programming 
(optimization) model developed for the force structure 
problem. Chapter IV presents the implementation of this model 
using a commercially available Linear Programming solver and 
a research level front-end user interface. Chapter V provides 
a sample problem using Fiscal Year 1989 data with analysis of 
the model's results and the impact of modifications to the 
various parameters used to describe the platforms. Chapter VI 
summarizes the findings of this study and proposes topics for 
further study. 
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II. MISSION REQUIREMENTS AND PLATFORM UTILIZATION 



As explained in Section B of Chapter I, previous Coast 
Guard projects or studies do not have a well defined 
methodology for quantifying mission requirements nor 
utilization of platforms. In this chapter, methods for 
quantifying these two factors in planning Coast Guard force 
structure are proposed. The main objective is to provide a 
reproducible and objective method for quantifying or 
estimating annual requirements for the seven major operating 
programs and specifying the usage of various platforms. Here, 
reproducibility means that, given the same set of data and 
using the proposed methods, different decision makers should 
be able to produce the same annual requirements and platform 
utilization. 

In addition, to gain acceptance in the Coast Guard 
community, new definitions or data requirements are kept to a 
minimum. Many parameters are provided to allow incorporation 
of judgmental factors which are best left to decision makers. 
In fact, officials at the Coast Guard headquarters in 
Washington, D.C. are kept abreast and express no objection to 
the proposed methodology. 
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A. QUANTIFYING ANNUAL MISSION REQUIREMENTS 

In the same major operating program, missions are 
generally different depending on the distance which must be 
traveled by the platform types, the length of the mission, the 
weather conditions, etc. Thus, it would be insufficient to 
simply specify total annual hours of platform usage for each 
operating program. When considering the fact that one 
platform may be capable to perform several types of missions, 
this insufficiency is more evident. For example, a platform 
designed primarily for ELT may not be able to perform a SAR 
mission at the same level of effectiveness and efficiency as 
a platform which is primarily designed for SAR. So, only 
specifying annual mission hours and allowing for multiple 
mission capability implies perfect substitution between 
platforms primarily designed for different types of missions. 

To allow for the difference in missions in the same 
program and more realistic degree of substitution, missions 
are furthered classified into three different classes. The 
choice and definition of the three classes are based on ten 
years of personal experience in the Coast Guard. However, 
they also serve to illustrate the basic concept. Decision 
makers with more experience and knowledge are more suitable to 
make the final determination. The implementation in Chapters 
IV and V also allow users to define their own classifications. 
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Table 2 and Table 3 display the classification of mission 
classes for surface and air platforms respectively. Among the 
three classes. Class I is the easiest and Class III is the 
hardest. Thus, missions in Class III should be mainly 
assigned to platforms specifically designed for the more 
difficult operating environment. On the other hand, missions 
in Class I may be performed by most if not all platforms in 
the Coast Guard inventory. The numerical attributes 
(Distance, Endurance, Speed, and Weather) are based on 
personal experience and they should be fully agreed upon by 
decision makers before full implementation. 



TABLE 2. 

MISSION CLASSES FOR SURFACE PLATFORMS 





Class I 


Class II 


Class III 


Distance 
off shore 


less than 
lOOnm 


100 - 200 nm 


200+ nm 


Endurance 


1 week 


1-2 weeks 


2+ weeks 


Speed 


at least lOkts 


lOkts - 20kts 


20+ kts 


Weather 


Good 


Fair 


Poor 



TABLE 3. 

MISSION CLASSES FOR AIR PLATFORMS 





Class I 


Class II 


Class III 


Distance 
off shore 


less than 
150nm 


150 - 300nm 


300 - 400nm 


Endurance 


3 hrs 


4 hrs 


5 hrs 


Weather 


Good 


Fair 


Poor 
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With this classification, the mission hours required by 
each major program are further subdivided into three classes, 
e.g., Class I-SAR, Class II-SAR and Class III-SAR. The next 
section then describes how to allocate time of each platform 
to each major operating program and class. 

B. PLATFORM UTILIZATION 

Although platforms may be multi-mission capable, they are 
generally designed to serve a specific program in a particular 
operating environment. Thus, care must be taken to prevent 
the over use of platforms in programs and operating 
environments not intended for them. Coast Guard documents 
such as the Sponsor's Requirements Document (SRD) and 
Commandant Instruction 3501.26 (series) serve as a foundation 
in establishing limits for the utilization of platforms to 
perform various missions in the three different classes. 

The Sponsor's Requirements Document details capabilities 
and limitations of platforms proposed for construction. 
Existing and new platforms must all have a SRD prepared prior 
to construction. It is based on this SRD that one can specify 
the maximum percent utilization of a given platform type in 
the three different classes of missions. Table 4 provides 
examples for maximum percent utilization for several platform 
types. 
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TABLE 4. 

CLASS UTILIZATION LIMITS 



PLATFORM 


CLASS I 


CLASS II 


- 

CLASS III 


82' WPB 


100% 


50% 


10% 


180' WLB 


100% 


100% 


50% 


378' WHEC 


100% 


100% 


100% 1 



The first row of Table 4 gives utilization limits for 82' 
WPB which is a patrol boat. From the limits, it can be 
ascertained that this patrol boat is primarily designed for 
mission of Class I type for it is allowed to use all (100%) of 
its available operating hours. However, it can only use 50% 
and 10% of its hours to perform missions in Class II and III 
respectively. Note that these percentages represent the 
maximum utilization in the three different classes; so they do 
not have to sum to 100%. 

Commandant Instruction 3501.26 (series) provides data 
concerning the Required Operational Capabilities (ROC) of 
various platforms. The format of this data is a table which 
lists the primary and secondary major operating programs 
suitable for each platform type. However, reviewing 
historical Abstract of Operations data reveals that there were 
platforms which performed missions not in their primary or 
secondary programs. This prompted an extension of the ROC 
categories for assigning maximum percent utilization of 
platforms to various programs. The basic idea is to allow a 
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high degree of utilization in the primary and secondary 
programs of each platform type. However, limited utilization 
would be allowed for missions not included in the primary or 
secondary programs. Table 5 gives examples of maximum 
percentage assignment for several platforms. 



TABLE 5. 

PROGRAM UTILIZATION LIMITS 



Platform 


SAR 


ELT 


ATON 


IOP 


MEP 


82' WPB 


80% 


100% 


10% 


0% 


0% 


180' WLB 


40% 


50% 


100% 


50% 


10% 


378' WHEC 


80% 


100% 


0% 


0% 


20% 



It can be gathered from the first row that ELT and SAR are 
the primary and secondary missions for an 82' patrol boat. 
The remaining 10% limit for ATON represent the utilization 
found in the historical Abstract of Operations data. Other 
programs with zero percent utilization limits are simply not 
suitable for an 82' WPB. 

Multiplying the class and program utilization tables 
together would provide the maximum percent utilization of a 
platform type by program and class. Although not explicitly 
stated, it should be clear that, prior to full implementation 
of the method described in this thesis, the percentages should 
be approved by all concerned parties in the Coast Guard. 
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Percentages used in this thesis are based on personal 
experience and offered only as an illustration. 

Chapter III describes the formulation of this problem as 
a mixed integer linear programming model. 
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III. MATHEMATICAL MODEL 



This chapter describes a mathematical programming model 
which allocates platforms to districts and local area commands 
to satisfy the annual mission requirements as quantified in 
Chapter II. The model attempts to use currently available 
platforms to satisfy the requirements. When necessary, the 
relocation of platforms from one district to another is 
permitted. However, the operational organization precludes 
the relocation of area platforms to districts. In addition, 
the model is also permitted to acquire new platforms if the 
available platforms are insufficient. It should be pointed out 
that it is not intended for the Coast Guard to purchase new 
platforms when the model's answer so states. Instead, 
allowing for new platforms is simply a tool to prevent model 
infeasibility and points out deficiency in the current 
platform availability. The decision to acquire new platforms 
requires models which allow for, e.g., multiple periods and 
service life extension programs. Models of this type are 
generally large and require data not yet available. Hence, 
consideration of such models is recommended for future 
investigation. 
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A. LINEAR MIXED INTEGER PROGRAM MODEL FORMULATION 



Indices: 




P 


platforms 


d, 6 = 


1,...,D districts 


m = 


1, . . . ,M missions 


1 


1 , . . . , 3 class 



Index Sets: 



G 


district surface platforms G c {1,...,P} 


$ 


district air platforms $ c {1,...,P} 


r 


area surface platforms T c (1,...,P) 

where {1,2,... ,P} = G u u T and Q, $ and T are 

mutually disjoint. 



Given and Derived Data: 



» P ..l 


upper utilization limits for hours expended on 
mission m at class 1 by platform p (see Chapter 
II, Section B) 


hr p 


maximum number of hours platform p can perform per 
year 


c p 


annual operating cost for platform p 


sh d. m .i 


annual surface hours required in district d for 
mission m at class 1 (see Chapter II, Section A) 




annual air hours required in district d for 
mission m at class 1 (see Chapter II, Section A) 


S Pm 


maximum percentage of air hours required which 
can be performed by surface platforms 


a Pm 


maximum percentage of surface hours required which 
can be performed by air platforms 


EP d 


current number of platform p located in district d 
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F p m ( factor for converting air hours to surface 
hours for mission m class 1 and by air 
platform p where p € $ 

G p m { factor for converting surface hours to air 
hours for mission m class 1 and by surface 
platform p where p € 0 u T 

a p cost of acquiring a new platform type p 

rP penalty for relocating platform p 

Nonnegative Variables: 

ZA p d m t number of surface hours performed by platform p in 
district d for mission m at class 1, where p e $ 
(see Figure 2) 

zs p d m ( number of air hours performed by platform p in 

district d for mission m at class 1, where p e 0 u 
T (see Figure 2) 

Y p d m ( number of hours platform p in district d performs 
mission m at class 1 (see Figure 2) 

A p d 8 number of district platforms p moved from district 
d to district 6 where p € Q u $ (see Figure 3) 




Figure 2. Representation of Hour Variables 
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Nonnegative Integer Variables: 



OXA p number of currently available (old) area 

surface platforms assigned where p e T 

OXD^ number of currently available (old) district 

platforms p assigned to district d 
where p € Q u # 

NXA P number of new area surface platforms assigned 

though not currently in inventory where p € T 

NX^ number of new district platforms p assigned to 

district d but not currently in inventory, 
where p e Q u 
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Formulation: 



MINIMIZE Y, E C*{OXD$ + NXD$) + Y C*(OXA p + NXA P ) 
d=l p€L?U<fc per 



+ E E aP wxDd + E 

d*i peUU 4 > per 



pel 
a p NXA p 



D D 



+ E E E rP 

d«i 6*1 pe^UQ 



Subject to: 




E Ff.i ZAf.m.i * sh d m l 

p€ 4> 



V d,m,l 



( 1 ) 



p€ 4> 



IdU.j + 



E ZS d,m.l * a ^d,m,l 

peiSur 



v a,/n. 



(2) 



E E Z ALm * a Pm E Sh d, ro ,i 

1-1 pea> J-l 



V ' d,m 



(3) 



E E z 5 d P -«.i * S P»E ah d.B.i V d,m 

l-l peOllr 1*1 



(4) 



M 3 



W 3 



E E + E E Z5 d p ^i * hi *{0X0$ + WXD d p ) V d,p e Q (5) 



JJl-l 1-1 D>-1 J-l 
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( 6 ) 



£ £ £ + £ £ E Z5 £».i * hr^OXA^ * NXAd) Vper 

d=l m=l 2-1 



DM 3 

EEL 

d-l /n-1 2-1 



M 3 M 3 

E £ *£«.i + £ £ * hr^OXDS + WXD d p ) V d,p e $> (7) 

in-1 2-1 jd-1 2-1 

Yd.m.i + ZSi, mtl <: hr pfpXDd + nxd£)u£j V d,m,l,p e fi (8) 



y£».j + 2AAJ..J * hr^OAPj’ + NXD§)u* J V d,m,l,p e <I> (9) 

r> 2 ? 

E *A».i + E z 5 d P ffi .i * hrP{OXDi * NXDgju*'! V m,l, P eT (10) 

d-l d-l 

g AS.* £ i? d P V d, p e QU$> (ll) 

]£A p dib = OXDf V 8,pefiU$ (12) 

d-l 

D 

OXA p ^Yj E d v per ( 13 ) 

d=i 



The objective function consists of three terms. The first 
term computes the annual operating cost of the platforms 
utilized by the model. The second and third terms calculate 
the acquisition and relocation cost, respectively. 
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Constraints (1) and (2) simply ensure that the number of 
hours performed by the various platforms meet the annual 
surface and air mission requirements, respectively. It should 
be noted that both constraints allow surface platforms to 
perform missions normally assigned to air platforms and vice 
versa. However, this cross utilization of platforms is not 
allowed to exceed certain percentages specified by the user. 
Constraint (3) then limits the number surface mission hours to 
be performed by air platforms to be no more than a percentage 
(ap m ) of the required surface mission hours. Similarly, 
constraint (4) places a limit on surface platforms. 

Constraints (5) through (7) ensure that there are enough 
platforms in each district and/or area to fulfill the hours 
allocated to perform various mission requirements as specified 
by constraints (1) through (4) . Constraints (8) through (10) 
guarantee that the upper utilization limits, as defined by 
Tables 4 and 5 in Chapter II, are not violated. These limits 
restrict the maximum number of hours that district and area 
platforms can perform missions in the three classes. 

Constraints (11) and (12) together account for the , 
relocation of existing platforms. Constraint (11) ensures 
that the number of platforms relocated from district d to 
other districts does not exceed the number of platforms 
originally stationed at d. Constraint (12) then guarantees 
that the number of platforms at district 6, relocated from 
other districts or otherwise, satisfies the allocation, 0XD 4 p , 
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specified in constraints (5) - (8) . Constraint (13) is 
included to ensure that the number of actual area platforms 
assigned does not exceed the current available number of area 
platforms. All platforms which exceed the current force 
structure must be considered new and be charged with the 
acquisition cost. 

The above model is basic to the force structure problem. 
A variety of embellishments are possible when additional data 
becomes available. As an example, constraints (14) and (15) 
below can be appended to the basic model to guarantee that the 
budgets for acquisition of new platforms and relocation of 
existing ones are not exceeded. 



D 




a p 



NXDg + £ a p NXA p £ alim 
per 



(14) 



£ E E rP * mlim 

d* l o^i pe DUO 



(15) 



In this case, the variables alim and mlim represent the two 
budgets. 

Chapter IV describes the implementation of this model and 
the development of a research level front-end user interface 
to facilitate in solving the force structure problem using the 
model developed above. 
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IV. IMPLEMENTATION AND USER INTERFACE 



The primary considerations for the selection of a LP 
solver and the programming language for the interface were 
portability and the current availability of software within 
the Coast Guard. The mathematical model described in Chapter 
III was implemented using the General Algebraic Modeling 
System (GAMS) [Ref. 13] available on the Amdahl 5990-500 
computer at the Naval Postgraduate School. To facilitate the 
input and modification of data for the model, a research level 
front-end user interface was written in FORTRAN. This 
interface allows for ease of data input and shields users away 
from the mathematical model implemented in GAMS. As a 
byproduct, the interface also helps in maintaining the data 
since it has many of the basic functions of a database system. 
The description of this interface is fully described below. 

A. USER INTERFACE 

To allow ease of use, an interface was developed to make 
the GAMS model transparent to the user. It is assumed that 
the user of this interface would be knowledgeable about the 
force structure problem but would not need to know the 
intricacies of the GAMS formulation or optimization 
techniques. The data required for this model is primarily 
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constant and should be store in data files to avoid needless 
repetitive data entry by the user. Storing and retrieving the 
model's data via files reduces the chance of data entry errors 
thus increasing the efficiency of the system and ultimately 
enhancing user acceptance. However, all data should be 
readily available to the user for modification as needed. 
Figure 4 provides a overview of the relationship between the 
data, the interface and GAMS. 




Figure 4: Data Flow Between User, Interface, and GAMS 



The interface is menu driven with a main menu which 
provides the user with access to all the functions shown in 
Figure 5. Most of these functional subroutines provide the 
user with different menus to further define the desire of the 
user. 
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Figure 5: Functional Routines of User Interface 

The user interface is a simple database program which 
allows modification of the data base, additions and deletions 
and the like. Additionally, the user must be able to view the 
database and produce output in either an easily understood 
format or in the format required by GAMS to solve the model. 
Though these two formats are not considered mutually 
exclusive, the ability to provide concise summary tables of 
the data for the user is preferable over printing out the 
sets, parameters, and tables used by the GAMS formulation. 
Figures 6 and 7 provide samples of the formatted output 
available from the interface's PRINT DATABASE option. 
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PLATFORM 


NAME: WPB110 




PLATFORM 


TYPE: DISTRICT SURFACE 






YEARLY OPERATING DATA: 




COST(SM): 0.590 HOURS: 


2000 


MISSION 


PROFILE DATA 




MISSION 


DESCRIPTION PERCENTAGE 


SAR 


SECONDARY MISSION 


50% 


ELT 


PRIMARY MISSION 


100% 


ATN 


NOT ASSIGNED 


0% 


IOP 


NOT ASSIGNED 


0% 


MEP 


TERTIARY MISSION 


10% 


MSAFE 


TERTIARY MISSION 


10% 


MILOP 


SECONDARY MISSION 


50% 


CLASS PROFILE DATA 




CLASS 


DESCRIPTION PERCENTAGE 


1 


FULLY CAPABLE 


100% 


2 


PARTIALLY CAPABLE 


50% 


3 


MARGINALLY CAPABLE 


10% 


CURRENT 


LOCATION OF ALL WPB110 PLATFORMS: 


FIRST 


2 FIFTH 2 SEVENTH 


16 


EIGHTH 


2 NINTH 0 





Figure 6: Sample Platform Database Output 
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DISTRICT: FIRST 

SURFACE HOURS REQUIRED AIR HOURS REQUIRED 
CLASS CLASS 



MISSION 


1 


2 


3 


1 


2 


3 


TOTAL 


SAR 


5000 


1000 


65 


668 


800 


400 


7933 


ELT 


8000 


6000 


1416 


800 


900 


2336 


19452 


ATN 


6000 


3000 


125 


200 


44 


0 


9369 


I0P 


800 


434 


200 


0 


0 


4 


1438 


HEP 


500 


48 


0 


328 


200 


0 


1076 


MSAFE 


500 


50 


0 


15 


0 


0 


565 


MILOP 


1500 


1000 


250 


800 


200 


133 


3883 


TOTAL 


22300 


11532 


2056 


2811 


2144 


2873 


43716 


CURRENT 


PLATFORMS LOCATED WITHIN 


FIRST 


DISTRICT BOUNDARIES: 


NAME 


QTY 


NAME 


QTY 


NAME 


QTY 


NAME 


QTY 



HU25 5 HH65 5 HH60 4 Cl 30 0 



WPB82 8 


WPB110 


2 


WTGB 


3 


WLB180 


3 


WLH157 2 


WLH133 


3 


UHEC 


1 


WMEC270 


2 



WMEC210 1 



Figure 7: Sample District Database Output 



The output is designed to provide the necessary 
information to the various decision makers involved with the 
Coast Guard force structure problem. The operational 
commanders would be most interested in the summary data 
concerning their district, the program managers would likewise 
be concerned with the data pertaining to the individual 
missions. The data displayed in Figures 6 and 7 must be 
formatted slightly differently to meet the structure of sets, 
tables and parameters required by GAMS. Figure 8 provides an 
example of the interface's ability to convert the database 
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into the data structures needed by GAMS for proper model 
definition and formulation. 



SETS 

P ALL PLATFORMS 

/ WPB82 , WPBl 10 ,WTGB / 

IDS(P) DISTRICT SURFACE PLATFORMS 
/ WPB82 , WPBl 10, UTGB / 

D DISTRICT 

/ FIRST , FIFTH , SEVENTH , EIGHTH , NINTH / 
M MISSION 

/ SAR ,ELT ,ATN ,IOP ,MEP ,MSAFE ,MILOP/ 



TABLE 

ADS(IDS,M) MAXIMUM PERCENTAGE OF TIME THAT PLATFORM (IDS) 

CAN 



PERFORM MISSION (M). 

SAR ELT ATN IOP 
WPB82 1.00 1.00 0.00 0.00 
WPBl 10 0.50 1.00 0.00 0.00 
UTGB 1.00 0.50 0.00 1.00 



MEP MSAFE MI LOP 
0.10 0.10 0.00 
0.10 0.10 0.50 
0.10 0.50 0.10; 



TABLE 

CURRENT(P,D) CURRENT FORCE 
FIRST FIFTH 
WPB82 8 

WPB110 3 2 

UTGB 3 1 



STRUCTURE FOR THE COAST GUARD 



SEVENTH EIGHTH NINTH 
9 10 0 

17 2 0 

0 5; 



PARAMETERS 

RHDS(IDS) MAXIMUM NUMBER OF RESOURCE HOURS PLATFORM (IDS) 
CAN PERFORM PER YEAR. 

/ WPB82 1650 

WPBl 10 2000 

UTGB 2300 / 

COSTDS(IDS) OPERATING COSTS FOR PLATFORM (IDS) PER YEAR IN 
MILLIONS OF DOLLARS. 

/ WPB82 0.370 

WPBl 10 0.590 

WTGB 0.650 /; 



Figure 8: Interface Output Following GAMS Data Structure 



Though the GAMS model representation is purported to be 
easily read by people and by computers [Ref. 13], the 
alternate output available through the interface provides the 
necessary information in a more compact and understandable 
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format for the user concerned more with the problem than with 
the formulation of the model. 

To prevent users from accidentally providing data for 
which no feasible solution exists, the interface performs a 
simple preprocessing of the entered data. This preprocessing 
checks that there is at least one platform type capable of 
performing each mission-class pair for which there is a 
positive mission hour requirement. Other forms of 
infeasibility would be detected by GAMS. 

B. GENERAL ALGEBRAIC MODELING SYSTEM 

The Coast Guard has recently acquired GAMS software for 
ongoing research into the force structure problem. The 
decision to implement this model using GAMS was motivated by 
the designed portability of the software [Ref. 13 :p. 4] and 
the Coast Guard's recent acquisition of this particular 
system. Another principle of GAMS design is that "the 
optimization problem should be expressible independently of 
the data it uses" [Ref. 13:p. 4]. Thus GAMS allows for the 
separation of the model's data from the actual formulation of 
the model. The data and formulation can be physically 
separated with the use of the INCLUDE statement. This feature 
allows the data for the model to be stored in a different file 
for later use with a GAMS formulation stored in another file. 
This also allows for numerous data files to be developed and 
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saved for use with the model while retaining only one copy of 
the model formulation and associated output statements. 

The FORTRAN-based solver called XMP/ZOOM [Ref. 13 & 14] is 
the only mixed integer programming solver currently available 
with the Coast Guard version of GAMS. By necessity, XMP/ZOOM 
was used as the integer programming solver for this 
mathematical model. During the development and testing of the 
mathematical model several problems became apparent with the 
use of XMP/ZOOM which is partly due to the fact that ZOOM "is 
intended for medium-sized problems with no special structure 
and up to about 200 zero/one variables" [Ref. 13:p. 225]. 

Chapter V will discuss the computational difficulties 
experienced with XMP/ZOOM and examine the feasibility of the 
optimization model by using Fiscal Year 1989 data mission hour 
requirements and platform availability. 
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V. SAMPLE PROBLEM 



After verifying the correctness of the model with a small 
data set, a realistic data set was considered in order to 
examine the operating characteristics of the model and its 
interface for the anticipated utilization by the Coast Guard. 
The sections below contain a description of how the data set 
was obtained, a discussion concerning the performance of the 
model and its solver, and sample analyses of the model 
outputs . 

A. PROBLEM DATA 

As stated in Chapter I, the Coast Guard organization is 
divided into two area commands, Atlantic and Pacific. Since 
they operate independently in that they do not share 
platforms, only the force structure problem for the Atlantic 
area command is considered. However, the Pacific area command 
can be modeled in a similar manner. The Atlantic area command 
contains six districts: First, Second, Fifth, Seventh, Eighth, 
and Ninth. The Second district consists of states in the 
midwest and is not adjacent to any coastal waters. This 
district employs only one type of platform, river buoy tenders 
to perform virtually one type of mission, ATON on the 
riverways. So, the problem of determining the number of 
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platforms for the Second district is simple and can be 
considered independently from the other districts in the 
Atlantic area command. Thus, the model contains one area 
command and five districts. 

Since this thesis only considers the force structure at 
the district and area levels, smaller platforms (vessels less 
than 65 feet in length) are excluded from the data set. Also 
excluded are icebreakers, construction tenders and the 
training vessel Barque Eagle. These platforms are used 
primarily for single missions and can be considered 
independently. Thus, the platform types included in this 
study are cutters and aircraft noted for their multi-mission 
design. The operating cost for these platforms were estimated 
from three sources: FY89 Standard Personnel Cost tables 
prepared by the Coast Guard Budget Division, cost tables 
provided by the Aviation Branch of the Coast Guard Operations 
Division, and discussions with Coast Guard officials. 
Additionally, to define utilization limits as described in 
Chapter II, the capabilities, design specifications, and 
manning levels for cutters and aircraft were extracted from 
Commandant Instruction 3501.26, The Ships and Aircraft of the 
U.S. Fleet , and personal experience. 

The required mission hours were extracted from the FY1989 
Abstract of Operations [Ref. 12 & 15]. Using information 
inferred from this document and personal experience, these 
mission hours were separated into the three mission classes as 
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defined in Chapter II. Mission hours in this data set include 
those hours utilized for missions benefiting the seven major 
operating programs; these hours do not include those used for 
training. To account for the required training hours, the 
available operating hours of each platform type are reduced by 
the number of hours needed for training (see Ref. 11) . 
Training is required for every platform and consumes a fairly 
constant amount of time. Therefore, subtracting the training 
hours from the available hours of each platform is not only a 
simple solution, but also prevents the model from assigning 
all training hours to, e.g., a single platform. 

Finally, acquisition costs and movement costs for the 
various platform types were set to constants, since these 
costs are not currently available. It should be emphasized 
that the data used in this sample problem is for illustrative 
purposes only. For example, required mission hours extracted 
from the Abstract of Operations would not be appropriate for 
actual planning, since these reports only document the mission 
hours performed by various platforms. In general, the total 
hours performed are restricted by the availability of 
platforms. When platforms are unavailable, some required 
missions are left unfulfilled and therefore not recorded in 
the Abstract of Operations . Thus the hours reported in this 
document underestimate the actual mission hours needed. 
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B. SOLVING THE MODEL 



During the development and testing of this model it was 
found that XMP/ZOOM consumed vast amounts of computer resource 
time while attempting to optimally solve the model. Given the 
uncertainty involved in determining future mission needs and 
the inaccuracy in forecasting operating costs for newly 
constructed platforms [Ref. 16], a near optimal solution, 
e.g., within twenty percent, would be acceptable. However, a 
ten percent or less tolerance is desirable. 

To further reduce the computational time to an acceptable 
level, the following cutting plane is added to the model. 

£ y) (OXD$ + NXD£) + Y^(OXA p + NXA p ) * T (16) 

d= l p€ lTG$ pit 

where T is the floor of the total number of platforms obtained 
by solving the model without the integrality restriction (or 
the relaxed model) . Theoretically, this cutting plane is not 
valid. However, based on the results in Table 6, the model 
with this cut closely approximated the original model. 

TABLE 6 

SOLUTIONS TO RELAXED MODELS WITH AND WITHOUT THE CUT 



Model 

Version 


Model w/o 
the Cut 


Model with 
the Cut 


Relative 

Difference 


I 


165.5081 


165.9127 


0.24% 


II 


158.6821 


158.7665 


0.05% 


III 


163.4569 


163.5397 


0.05% 


IV _ 


169.1057 


169.1797 


0.04% 
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Table 6 reports the objective function values for four 
versions of the relaxed model with and without the cut. Since 
the two sets of values do not differ significantly it is 
assumed that the feasible region of the integer model with and 
without the cut are equivalent or nearly so. 

Without the cut, version II of the model consumed over 
five hours of CPU time and still failed to reach an integer 
solution. With the cut, the model used 30 minutes of CPU time 
to produce an integer solution within 15 percent of 
optimality. On the average, models with the cut took 
approximately 23 minutes to obtain an acceptable (within 15 
percent) integer solution. Thus, to perform sample analyses, 
the cut was added to each version of the model discussed 
below. 

C. OUTPUTS AND ANALYSIS 

Figure 9 displays the current locations of all district 
and area platforms used for all runs mentioned in this 
section. These locations were obtained, as before, from 
FY1989 Abstract of Operations. To insure that data in Figure 
9 and those gathered in Section B yield a meaningful problem, 
the model in Chapter III was solved with the cross utilization 
factors of platforms set to zero. This disallows the surface 
platforms to perform air missions and vice versa. The 
resulting model contains 2437 continuous variables, 106 
discrete variables, and 1182 equations and took 1,758 CPU 
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seconds to obtain a solution within seven percent of 
optimality. The result from the ZOOM solver showed that the 
data set is meaningful for it generates a feasible problem. 



PLATFORM 


1ST 


5TH 


DISTRICT 
7TH 8TH 


9TH 


AREA 


TOTAL 


AIRCRAFT 


HU 2 5 


5 


0 


10 


3 


0 




18 


HH65 


5 


3 


16 


11 


4 


- 


39 


HH3 


4 


3 


9 


0 


2 


- 


18 


C130 


0 


5 


8 


0 


0 




13 


SURFACE 


WPB82 


8 


8 


8 


10 


0 




34 


WPB110 


2 


2 


16 


2 


0 




22 


WTGB 


3 


1 


0 


0 


5 


- 


9 


WLB180 


3 


3 


2 


2 


5 




15 


WLM157 


2 


3 


0 


0 


0 




5 


WLM133 


3 


0 


1 


2 


0 




6 


WHEC 


- 


- 


- 


- 


- 


1 


1 


WMEC270 


- 


- 


- 


- 


- 


10 


10 


WMEC210 












10 


10 



Figure 9: Fiscal Year 1989 Coast Guard Force Structure 



Figures 10 - 12 are part of the outputs from the model. 
A full listing of model outputs are provided in Appendix C. 
Figure 10 shows how the model reallocates the platforms. 
Recall that area platforms can not be relocated to districts 
and vice versa. So, area platforms are not relocated and 
district platforms must be relocated among themselves. It is 
interesting to note that the model does not utilize all 
available platforms. Figure 11 displays the model allocation 
of the operating hours of each platform type to each major 
operating program. To provide a more detail allocation of 
these operating hours, the model also provides the allocation 
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by districts. For example, Figure 12 illustrates the 
allocation for the seventh district. 



PLATFORM 


DISTRICT 












1ST 


5TH 


7TH 8TH 


9TH 


AREA 


TOTAL 


AIRCRAFT 
















HU25 


1 


1 


2 0 


0 


- 


4 




HH65 


7 


6 


16 8 


2 


- 


39 




HH3 


4 


0 


12 0 


0 


- 


16 




C130 


0 


1 


2 0 


1 


— 


4 




SURFACE 
















WPB82 


8 


6 


15 5 


0 


- 


33 




WPB110 


2 


2 


16 2 


0 


- 


22 




WTGB 


3 


1 


1 0 


4 


- 


9 




WLB180 


0 


0 


1 2 


2 


- 


5 




WLM157 


2 


3 


0 0 


0 


- 


5 




WLM133 


1 


1 


2 2 


0 


- 


6 




WHEC 


- 






- 


0 


0 




WMEC270 


- 


- 


- - 


- 


5 


5 




WMEC210 


— 


— 


— “ 


— 


10 


10 




Figure 10: 


Model ' 


s Proposed 1989 Coast Guard Force Structure 


PLATFORM 






MISSION 










SAR 


ELT 


ATN 


IOP 


MI LOP 


MSAFE 


MEP 


AIRCRAFT 
















HU25 


2960 


1427 


123 


4 


1020 


15 


50 


HH65 


2158 


17069 


603 


74 


418 


47 


998 


HH3 


4250 


4891 






183 


16 


170 


C130 


1723 


2048 






229 






SURFACE 
















WPB82 


11382 


32073 








738 


357 


WPB110 


993 


36762 






6195 




50 


WTGB 


2875 


2233 




2992 


415 


136 




WLB180 


681 


2121 


23269 


1800 


123 




5 


WLM157 




491 


11508 










WLM133 




1001 


16450 


472 


770 




548 


WMEC270 


410 


34150 












WMEC210 




42558 






641 







Figure 11: Model Allocation of Platform Hours 
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PLATFORM 






MISSION 








SAR 


ELT 


ATN 


MILOP 


MSAFE 


MEP 


AIRCRAFT 














HU25 


2221 


900 


79 








HH65 




10320 










HH3 


2296 


2800 






16 


170 


C130 


123 


2048 




229 






SURFACE 














WPB82 


5200 


7848 






66 


86 


WPB110 


255 


30262 




1483 






WLB180 




669 


7331 








WLM133 




231 


3619 








WMEC270 




27171.5 










WMEC210 




38868.5 











Figure 12: Allocation of Platform Hours for Seventh District 



Given that the above data set is meaningful, two types of 
analysis of interest to the Coast Guard were performed. The 
first one concerns how a change in the difficulty of the 
mission requirements effect the reallocation of platforms. 
The other concerns a similar effect due to changes in 
requirements of a mission type. As before, it is also assumed 
that there is no cross utilization between air and surface 
platforms. This assumption is necessary because of the fact 
that the Abstract of Operations did not record the cross 
utilization of air and surface platforms and that ZOOM would 
take an inordinate amount of time otherwise. 

The assumption of no cross utilization allows the model to 
be decomposed into two smaller problems, aircraft allocation 
and surface platform allocation. The independently solved 
problems can be combined to achieve the complete solution. 
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The original data set was used once again to examine the 
effects of this decomposition. The sum of the resulting two 
sub-problems contained 2,721 continuous variables, 44 discrete 
variables, and 1,181 equations and took 72 CPU seconds to 
obtain a solution within eight percent of optimality. This 
eight percent optimality yields an optimal objective function 
value of $170.76 million instead of $169.07 million from the 
run with seven percent optimality. Figure 13 indicates the 
combined force structure proposed by the two sub models. This 
95% reduction in computational time provided the necessary 
speed to complete the following two analyses. 



PLATFORM 


1ST 


5TH 


DISTRICT 
7TH 8TH 


9TH 


AREA 


TOTAL 


AIRCRAFT 


HU25 


1 


0 


2 


0 


0 


- 


3 


HH65 


7 


6 


16 


8 


2 


- 


39 


HH3 


3 


1 


13 


0 


1 


- 


18 


C130 


1 


1 


1 


0 


0 


— 


3 


SURFACE 


WPB82 


8 


5 


15 


5 


0 


- 


33 


WPB110 


2 


2 


16 


2 


0 


- 


22 


WTGB 


3 


1 


1 


0 


4 


- 


9 


WLB180 


0 


2 


3 


0 


3 


- 


8 


WLM157 


2 


3 


0 


0 


0 




5 


WLM133 


1 


0 


1 


4 


0 


- 


6 


WHEC 


- 


- 


- 


- 


- 


0 


0 


WMEC270 


- 


- 


- 


- 


- 


4 


4 


WMEC210 




— 


— 


— 




10 


10 


Figure 13 : 


Models 


1 Proposed 


Total Force 


Structure. 




Chancrina 


the 


mission 


difficulty: 


To 


examine 


how the 



change in the difficulty of the missions effects the model 
solution, the required mission hours were modified as follows. 
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Ten percent of required hours from Class II mission were 
shifted to Class III. Similarly, ten percent from Class I 
were also shifted to Class II. Recall that Class I mission 
hours are the easiest and Class III hours the most difficult. 
Therefore, the effect of this shifting of hours was to make 
the overall hours more difficult to perform. Since some hours 
were removed from Class II and some hours were added to Class 
II, the total percentage of hours changed was 9.9 percent, 
slightly less than the ten percent shifted between the 
classes. Figure 14 provides the revised force structure 



proposed by the model. 



PLATFORM 








DISTRICT 








1ST 


5TH 


7TH 


8TH 


9TH 


AREA 


TOTAL 


AIRCRAFT 


HU25 


1 


0 


2 


0 


1 


- 


4 


HH65 


7 


6 


16 


8 


2 


- 


39 


HH3 


4 


1 


13 


0 


0 


- 


18 


C13 0 


0 


1 


1 


0 


0 


— 


2 


SURFACE 


WPB82 


8 


5 


13 


5 


0 


- 


31 


WPB110 


2 


2 


16 


2 


0 


- 


22 


WTGB 


3 


1 


0 


0 


5 


- 


9 


WLB180 


1 


1 


0 


1 


2 


- 


5 


WLM157 


0 


3 


1 


1 


0 


- 


5 


WLM133 


2 


0 


2 


2 


0 


- 


6 


WHEC 


- 


- 


- 


- 


- 


1 


1 


WMEC270 


- 


- 


- 


- 


- 


6 


6 


WMEC210 












10 
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Figure 14 : 


Force 


Structure For 


More 


Difficult 


Missions. 



Numerous changes to the force structure were evident and 
are summarized below in Figure 15. 
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Figure 15: 
Difficulty. 



Changes in Force Structure Due to Increased 



Interestingly, the overall number of platforms required to 
meet the new mission hour requirements actually decreased by 
two. However, the platforms not allocated were generally 
those designed for easier mission hours and consequently, less 
expensive than the more capable platforms. The number of area 
platforms, considered the most capable for a wide range of 
missions, increased by three due to the change in mission 
requirements. Considering that the force structure problem is 
concerned with the location of platforms as well as the total 
number of platforms, the shifting of platforms from one 
district to another must be considered a change in the force 
structure. By this definition of change, the 9.9 percent 
change in mission difficulty resulted in 26 changes or a 16.25 
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percent change in the force structure. This indicates small 
changes in difficulty of the missions produces significant 
changes in the force structure. This emphasizes the necessity 
of having accurate required mission hours. 

Increasing required number of ELT mission hours : To 

examine the impact of increased operations within a particular 
mission. ELT was chosen for this example since its required 
mission hours accounted for over half of the total mission 
hours accomplished by the Coast Guard in Fiscal Year 1989. 
Figures 16 and 17 graphically illustrate changes in the force 
structure caused by this increase in required ELT mission 
hours. 




Figure 16: Additional Platforms Grouped by Program Rating 
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Figure 17: Additional Platforms Allocated to Districts 



Recall when using the original data set of required 
mission hours, the model does not utilize all available 
platforms. As such, the additional platforms described along 
the y-axis of Figures 16 and 17 are not new acquisitions, 
rather they are platforms previously not utilized by the 
model. In this case, a 40 percent increase in required ELT 
mission hours or, equivalently, 27.2 percent increase in 
overall mission hours, resulted in only a 15 percent increase 
in the number of platforms. This increase in the number of 



46 



platforms was relatively small because of the multi-mission 
capability of the platforms. 

Figure 16 depicts the increase in the utilization of 
platforms separated into two groups. The first group consists 
of platforms primarily designed for ELT as indicated in ROC 
documentation. The other consists of the remaining platforms. 
Clearly, this figure shows a marked increase of approximately 
one platform for every two percent increase in mission hours 
for platforms primarily designed for ELT. However, the 
remaining platforms increased at a slower rate, approximately 
one for every six percent. One conclusion for this behavior 
is that platforms are now being used for their primary mission 
which in this case is ELT. 

Force structure changes in the Seventh district was of 
interest since originally 68 percent of the required ELT 
mission hours were in this district. Figure 17 compares the 
increase in district platforms allocated to both the Seventh 
district and the remaining Atlantic area districts. As 
expected, the number of platforms increased more rapidly in 
this district because it has the largest percentage of 
required ELT mission hours. At a 15 percent increase in 
required ELT hours, those area platforms previous unused by 
the model were being allocated to accomplish the increasing 
number of ELT mission hours. At this point two additional 
area platforms were allocated solely for accomplishing ELT 
hours in the Seventh district. After the 15 percent increase 
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one additional area platform was allocated for each five 
percent increase in required mission hours. 

Examining figures similar to Figure 16 and 17 can provide 
insight into the adequacy of the available platforms in 
response to an unexpected increase in mission hours such as 
the Cuban Boat Lift of 1980 or hurricane HUGO in 1990. Other 
analyses are possible by examining other portions of the 
model's output. This section illustrates only a few of the 
analytical possibilities available. 
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VI . CONCLUSION 



This thesis developed and implemented an optimization 
based tool to aid the Coast Guard in establishing its force 
structure. Since the Coast Guard has questioned previous 
methods of defining mission requirements and platform 
utilization. Chapter II proposed new methods to systematically 
quantify these factors. Then, Chapter III specified an 
integer programming model to allocate platforms to meet the 
requirements in the manner specified by the utilization 
criteria in Chapter II. 

A user friendly interface was developed in Chapter IV. 
This interface is designed to help users maintain and update 
data for the problem as well as to shield them from the 
mathematical programming aspects of the problem. However, the 
interface still provides sufficient flexibility for the user 
to explore different but related versions of the basic force 
structure problem. Finally, sample outputs and analyses using 
the implemented system were given in Chapter V. 

A. AREAS OF FUTURE STUDY 

This study provides a foundation with which the Coast 
Guard can develop a powerful tool to solve its force structure 
problem. Future investigations in this direction would not 
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only benefit the Coast Guard in increasing its efficiency and 
effectiveness, but are also necessary when facing an 
unfavorable budgetary environment. Among many future topics, 
three are presented here. 

(1) As discussed earlier accurate input data is necessary 
to improve the usefulness of this model. Toward this end, the 
Coast Guard is encouraged to investigate the development of 
models or improved techniques for forecasting costing data and 
mission requirements data. With improvements to the input 
data, the ability to develop a multi-period model becomes a 
logical extension to this study. A multi-period model would 
then be able to look at such issues as hull life, the 
congressionally mandated acquisition process as defined in OMB 
Circular Number A-109, and service life extension programs 
(SLEP) such as Fleet Rehabilitation and Maintenance (FRAM) and 
Mid-life Maintenance Availability (MMA) . All of these issues 
could then be used to determine the life cycle cost of each 
platform vice the simple annual costs currently used in this 
model, providing a more realistic view of the actual cost of 
a platform type. 

(2) Possible errors within the GAMS/ZOOM software were 
discovered during the implementation of the mathematical model 
and have been documented in previous studies using GAMS/ZOOM 
[Ref. 17]. The Coast Guard should investigate the possibility 
of procuring another commercial mixed integer solver for use 
with GAMS or the development of a customized problem generator 
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and integer solver written in a general-purpose compiled 
language (e.g., FORTRAN). 

(3) The research level front end interface is admittedly 
rudimentary and was developed to provide an example of the 
ability to insulate the user from an underlying mathematical 
model. An improved interface would greatly enhance the 
usability and acceptance of the optimization model. The Coast 
Guard is presently investigating a KSS (knowledge-based 
decision support systems) project using a DSS (decision 
support system) environment [Ref. 6]. The "DSS is an 
interactive software tool for working with models and data" 
[Ref. 6:p. 6]. The combination of this mathematical model 
with an interface based on the KSS environment would provide 
the Coast Guard with a necessary tool for analyzing the many 
force structure issues. 

In closing, it should be mentioned that this study 
provides the Coast Guard with one method of addressing its 
force structure problem. Since no one, mathematical or 
otherwise, model can completely capture every facet of the 
problem, the Coast Guard is also encouraged to consider other 
approaches, e.g., stochastic or statistical modeling, along 
with directions mentioned above. 
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APPENDIX A 



INTERFACE OUTPUT: GAMS DATA SET 



SETS 

P ALL PLATFORMS 

/ 

HU25 , HH65 ,HH60 ,C130 ,WPB82 ,WPB110 , 

WTGB , WLB180 ,WLM157 ,WLM133 ,WHEC ,WMEC270 , 

WMEC210 / 

IDS (P) DISTRICT SURFACE PLATFORMS 

/ 

WPB82 , 

WPB110 , 

WTGB , 

WLB180 , 

WLM157 
WLM133 / 

IDA(P) DISTRICT AIR PLATFORMS 

/ 

HU25 , 

HH65 , 

HH60 , 

C13 0 / 

IAS (P) AREA SURFACE PLATFORMS 

/ 

WHEC , 

WMEC270 , 

WMEC210 / 

D DISTRICT 

/ 

FIRST , FIFTH , SEVENTH , EIGHTH , 

NINTH / 

M MISSION 

/ SAR , ELT , ATN , IOP , MILOP, MSAFE , MER / 

L CLASS /l, 2, 3/; 

SET INDP(D) INDEPENDENT DISTRICTS / NINTH /? 
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ALIAS (D, DP) ; 
TABLE 





ADS ( IDS, M) 


MAXIMUM : 


PERCENTAGE OF TIME 


THAT PLATFORM 


* 




CAN PERFORM MISSION (M) . 










SAR ELT 


ATN 


IOP MILOF 


• MSAFE 


MER 




WPB82 


1.00 1.00 


0.00 


0.00 0.00 


0.10 


0.10 




WPB110 


0.50 1.00 


0.00 


0.00 0.50 


0.50 


0.10 




WTGB 


1.00 0.50 


0.00 


1.00 0.10 


0.50 


0.00 




WLB180 


0.50 0.10 


1.00 


0.10 0.10 


0.10 


0.10 




WLM157 


0.50 0.10 


1.00 


0.10 0.10 


0.00 


0.10 




WLM133 


0.50 0.10 


1.00 


0.10 0.10 


0.00 


0.10 


TABLE 














BDS (IDS , L) 


MAXIMUM ; 


PERCENTAGE OF TIME 


THAT PLATFORM 


* 




CAN PERFORM AT 


CLASS (L) . 










1 


2 


3 








WPB82 


1.00 


0.10 


0.00 








WPB110 


1.00 


0.50 


0.10 








WTGB 


1.00 


0.50 


0.10 








WLB180 


1.00 


1.00 


0.50 








WLM157 


1.00 


1.00 


0.10 








WLM133 


1.00 


0.50 


0.10 




• 

9 


TABLE 














ADA ( I DA , M ) 


MAXIMUM 


PERCENTAGE OF TIME 


THAT PLATFORM 


* 




CAN PERFORM MISSION (M) . 










SAR ELT 


ATN 


IOP MILOP MSAFE 


MER 




HU 2 5 


1.00 1.00 


0.10 


0.10 0.50 


0.10 


0.50 




HH65 


1.00 1.00 


0.50 


0.50 0.50 


0.10 


0.50 




HH60 


1.00 1.00 


0.50 


0.00 0.50 


0.10 


0.50 




C130 


1.00 1.00 


0.50 


0.50 0.50 


0.50 


0.50 


TABLE 














BDA (IDA, L) 


MAXIMUM 


PERCENTAGE OF TIME 


THAT PLATFORM 


* 




CAN PERFORM AT 


CLASS (L) . 










1 


2 


3 








HU25 


1.00 


1.00 


0.10 








HH65 


1.00 


0.50 


0.00 








HH60 


1.00 


1.00 


0.50 








C13 0 


1.00 


1.00 


1.00 




• 

9 


TABLE 














AAS ( IAS , M) 


MAXIMUM 


PERCENTAGE OF TIME 


THAT PLATFORM 


* 




CAN PERFORM MISSION (M) . 










SAR ELT 


ATN 


IOP MILOP MSAFE 


MER 




WHEC 


0.50 1.00 


0.00 


0.00 1.00 


0.00 


0.10 




WMEC270 


0.50 1.00 


0.00 


0.00 1.00 


0.50 


0.00 




WMEC210 


0.50 1.00 


0.00 


0.00 0.50 


0.50 


0.10 
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TABLE 



BAS ( IAS , L) 


MAXIMUM 


PERCENTAGE 


OF TIME THAT PLATFORM (IAS) 


* 


CAN PERFORM AT CLASS (L) . 




1 


2 


3 


WHEC 


1.00 


1.00 


1.00 


WMEC270 


1.00 


1.00 


0.50 


WMEC210 


1.00 


1.00 


0.50 

• 

9 



TABLE 



HRSURF ( D , 


M, L) SURFACE RESOURCE HOURS 


REQUIRED FOR MISSION 




(M) 


IN DISTRICT 


(D) AT 


CLASS (L) PER YEAR. 






1 


2 


3 


FIRST 


. SAR 


5000 


1000 


65 


FIRST 


.ELT 


8000 


6000 


1416 


FIRST 


. ATN 


5000 


3000 


1125 


FIRST 


• IOP 


800 


434 


200 


FIRST 


.MILOP 


1500 


1000 


250 


FIRST 


. MSAFE 


500 


50 


0 


FIRST 


.MER 


500 


48 


0 


FIFTH 


.SAR 


1000 


1200 


61 


FIFTH 


.ELT 


7000 


4800 


1045 


FIFTH 


.ATN 


6000 


4000 


1314 


FIFTH 


.IOP 


0 


0 


0 


FIFTH 


.MILOP 


1100 


1065 


900 


FIFTH 


. MSAFE 


59 


0 


0 


FIFTH 


.MER 


200 


71 


0 


SEVENTH 


.SAR 


4000 


1200 


255 


SEVENTH 


.ELT 


40000 


35050 


30000 


SEVENTH 


.ATN 


10000 


900 


50 


SEVENTH 


.IOP 


0 


0 


0 


SEVENTH 


.MILOP 


1000 


283 


200 


SEVENTH 


. MSAFE 


66 


0 


0 


SEVENTH 


.MER 


86 


0 


0 


EIGHTH 


.SAR 


1800 


620 


0 


EIGHTH 


.ELT 


6090 


7700 


4000 


EIGHTH 


.ATN 


12700 


1070 


0 


EIGHTH 


.IOP 


0 


0 


0 


EIGHTH 


.MILOP 


512 


150 


150 


EIGHTH 


.MSAFE 


63 


0 


0 


EIGHTH 


.MER 


50 


0 


0 


NINTH 


.SAR 


140 


0 


0 


NINTH 


.ELT 


289 


0 


0 


NINTH 


.ATN 


2300 


2000 


1769 


NINTH 


.IOP 


600 


2000 


1230 


NINTH 


.MILOP 


35 


0 


0 


NINTH 


.MSAFE 


136 


0 


0 


NINTH 


• MER 


5 


0 


0 
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TABLE 



HRAIR (D, M, 


L) AIR 


RESOURCE HOURS 


REQUIRED 


FOR MISSION 




IN 


DISTRICT (D) AT 


CLASS (L) 


PER YEAR. 






1 


2 


3 


FIRST 


• SAR 


668 


800 


400 


FIRST 


.ELT 


800 


1900 


1336 


FIRST 


. ATN 


200 


44 


0 


FIRST 


.IOP 


0 


0 


4 


FIRST 


.MILOP 


600 


300 


233 


FIRST 


. MSAFE 


10 


5 


0 


FIRST 


.MER 


278 


200 


50 


FIFTH 


.SAR 


900 


600 


360 


FIFTH 


.ELT 


1070 


700 


500 


FIFTH 


.ATN 


100 


71 


0 


FIFTH 


.IOP 


0 


1 


0 


FIFTH 


.MILOP 


121 


100 


70 


FIFTH 


.MSAFE 


18 


0 


0 


FIFTH 


.MER 


149 


100 


0 


SEVENTH 


.SAR 


2300 


1300 


1040 


SEVENTH 


.ELT 


5600 


5300 


5168 


SEVENTH 


.ATN 


79 


0 


0 


SEVENTH 


.IOP 


0 


0 


0 


SEVENTH 


.MILOP 


60 


69 


100 


SEVENTH 


.MSAFE 


10 


6 


0 


SEVENTH 


.MER 


120 


50 


0 


EIGHTH 


.SAR 


800 


870 


0 


EIGHTH 


.ELT 


1600 


1420 


0 


EIGHTH 


.ATN 


63 


0 


0 


EIGHTH 


.IOP 


0 


0 


0 


EIGHTH 


.MILOP 


100 


61 


0 


EIGHTH 


. MSAFE 


10 


4 


0 


EIGHTH 


.MER 


106 


100 


0 


NINTH 


.SAR 


700 


354 


0 


NINTH 


.ELT 


42 


0 


0 


NINTH 


• ATN 


150 


19 


0 


NINTH 


.IOP 


50 


23 


0 


NINTH 


.MILOP 


26 


10 


0 


NINTH 


.MSAFE 


12 


3 


0 


NINTH 


.MER 


55 


0 


0 



TABLE 



CURRENT (P,D) CURRENT FORCE STRUCTURE FOR THE COAST GUARD 





FIRST 


FIFTH 


SEVENTH 


EIGHTH 


NINTH 


HU25 


5 


0 


10 


3 


0 


HH65 


5 


3 


16 


11 


4 


HH60 


4 


3 


9 


0 


2 


C13 0 


0 


5 


8 


0 


0 


WPB82 


8 


8 


8 


10 


0 


WPB110 


2 


2 


16 


2 


0 


WTGB 


3 


1 


0 


0 


5 


WLB180 


3 


3 


2 


2 


5 


WLM157 


2 


3 


0 


0 


0 
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WLM133 

WHEC 

WMEC270 

WMEC210 



3 

1 

2 

1 



0 

0 

5 

3 



1 

0 

3 

4 



2 

0 

0 

2 



0 

0 

0 

0 



SCALARS 

PICK FLAG FOR WHICH OBJECTIVE FUNCTION TO USE 

CNT PLATFORM COUNT 

ACQLIM FLAG FOR LIMITING ACQUISITION COSTS 

MOVLIM FLAG FOR LIMITING MOVEMENT COSTS 

LIMACQ LIMIT FOR ACQUISITION COSTS 

LIMMOV LIMIT FOR MOVEMENT COSTS 

PENALTY 1 ACQUISITION COST CONSTANT 
PENALTY 2 MOVEMENT COST CONSTANT 

NOIDS FLAG FOR NO DISTRICT SURFACE PLATFORMS DEFINED 

NOIAS FLAG FOR NO AREA SURFACE PLATFORMS DEFINED 

NOIDA FLAG FOR NO DISTRICT AIR PLATFORMS DEFINED 



PARAMETERS 



AIRPCT(M) 



% OF SURFACE RESOURCE HOURS AIRCRAFT CAN 
ACCOMPLISH 



/ SAR 


0.00 




ELT 


0.00 




ATN 


0.00 




IOP 


0.00 




MILOP 


0.00 




MSAFE 


0.00 




MER 


0.00 


/ 


SURPCT(M) % 


OF AIR 


RESOURCE 


/ SAR 


0.00 




ELT 


0.00 




ATN 


0.00 




IOP 


0.00 




MILOP 


0.00 




MSAFE 


0.00 




MER 


0.00 


/ 


RHDS ( IDS ) MAXIMUM 


NUMBER OF 



CAN PERFORM PER YEAR. 



/ 



WPB82 


1650 


WPB110 


2000 


WTGB 


2300 


WLB180 


4000 


WLM157 


4000 
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WLM133 



3850 / 



RHDA(IDA) MAXIMUM NUMBER OF RESOURCE HOURS PLATFORM (IDA) 
CAN PERFORM PER YEAR. 

/ HU25 800 

HH65 645 

HH60 700 

C130 800 / 

RHAS (IAS) MAXIMUM NUMBER OF RESOURCE HOURS PLATFORM (IAS) 
CAN PERFORM PER YEAR. 



/ 


WHEC 


4320 






WMEC270 


4320 






WMEC210 


4320 


/ 


COSTDS (IDS ) 


OPERATING 


COSTS FOR 


PLATFORM 




MILLIONS 


OF DOLLARS. 




/ 


WPB82 


0.370 






WPB110 


0.590 






WTGB 


0.650 






WLB180 


1.580 






WLM157 


1.130 






WLM133 


0.990 


/ 


COSTDA (IDA) 


OPERATING 


COSTS FOR 


PLATFORM 




MILLIONS 


OF DOLLARS. 




/ 


HU2 5 


1.640 






HH65 


1.110 






HH60 


1.470 






C130 


2.220 


/ 


COSTAS (IAS) 


OPERATING 


COSTS FOR 


PLATFORM 




MILLIONS 


OF DOLLARS. 




/ 


WHEC 


4.280 






WMEC270 


2.690 






WMEC210 


2.130 


/ 



TOTCURRENT ( P) TOTAL AVAILABLE PLATFORMS 



ULDS ( I DS , M , L) 



ULAS ( IAS , M, L) 



ULDA(IDA,M, L) 



UPPER UTILIZATION LIMITS FOR RESOURCE HOURS 
EXPENDED ON MISSION (M) AT CLASS (L) 

FOR DISTRICT SURFACE PLATFORMS (IDS) 

UPPER UTILIZATION LIMITS FOR RESOURCE HOURS 
EXPENDED ON MISSION (M) AT CLASS (L) 

FOR AREA SURFACE PLATFORMS (IAS) 

UPPER UTILIZATION LIMITS FOR RESOURCE HOURS 
EXPENDED ON MISSION (M) AT CLASS (L) 

FOR DISTRICT AIR PLATFORMS (IDA) ; 



ULDS (IDS ,M,L) = ADS (IDS ,M) *BDS (IDS , L) ; 
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ULAS ( IAS , M , L) = AAS ( IAS , M) *BAS ( IAS , L) ; 
ULDA(IDA,M,L) = ADA(IDA,M) *BDA(IDA, L) ; 



RHDS (IDS ) = RHDS (IDS)/100.0; 
RHDA(IDA) = RHDA ( I DA )/100.0; 

RHAS (IAS) = RHAS (IAS)/ 100. 0; 

HRSURF ( D , M , L) = HRSURF (D,M, L) /100 . 0 ; 
HRAIR(D,M,L) = HRAIR (D,M, L) /100 . 0 ; 
TOTCURRENT ( P) = SUM(D, CURRENT (P, D) ) ; 

NOIDS = 1; 

NOIDA = 1; 

NOIAS = 1; 



PICK = 1.0; 
ACQLIM = 
MOVLIM = 
LIMACQ = 
LIMMOV = 
PENALTY 1 = 
PENALTY 2 = 



1 . 0000 ; 

1 . 0000 ; 

0 . 0000 ; 

0 . 0000 ; 

0.5000; 

o.iooo; 
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APPENDIX B 



GAMS FORMULATION 



$TITLE United States Coast Guard Fleet-Mix Problem 

$STITLE LT J.E. TOMKO , USCG 

* GAMS and Dollar Control Options 



$OFFUPPER OFFSYMXREF OFFSYMLIST 

OPTIONS LIMCOL = 0, LIMROW = 0, SOLPRINT = OFF; 
OPTIONS RESLIM = 50000, ITERLIM = 2500000; 
OPTIONS WORK = 35000; 

* Definitions and Data 



OPTIONS OPTCR = 0.20; 



$ INCLUDE GAMS DATA A 



* Model 

VARIABLES 

TCOST total operating cost of fleet per year; 

INTEGER VARIABLES 

NXDS(IDS,D) new # of district surface platforms (IDS) 

in district (D) 

NXAS(IAS) new # of area surface platforms (IAS) 

NXDA (IDA, D) new # of district air platforms (IDA) 

in district (D) 

OXDS (IDS , D) current # of district surface platforms (IDS) 

in district (D) 

OXAS (IAS) current # of area surface platforms (IAS) 

OXDA (IDA, D) current # of district air platforms (IDA) 

in district (D) ; 



OXDS. UP ( IDS, D) = TOTCURRENT (IDS) ; 

OXDA. UP ( IDA, D) = TOTCURRENT (IDA) ; 

OXAS. UP (I AS) = TOTCURRENT ( IAS ) ; 

POSITIVE VARIABLES 

YDS (IDS, D,M, L) # hrs of platform (IDS) in district (D) assigned 

to mission area (M) at class (L) performed 
for surface missions 

YAS ( IAS , D , M , L) # hrs of platform (IAS) in district (D) assigned 

to mission area (M) at class (L) performed 
for surface missions 
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YDA (IDA, D,M, L) 



MS (IDS , D, DP) 
MA(IDA, D, DP) 
ZDS (IDS , D, M, L) 



ZAS ( IAS , D , M , L) 



ZDA(IDA,D,M,L) 



# hrs of platform (IDA) in district (D) assigned 
to mission area (M) at classs (L) performed 
for air missions 

movement variable for platform type (IDS) 

movement variable for platform type (IDA) 

# hrs of platform (IDS) in district (D) assigned 
to mission area (M) at class (L) performed 
for air missions 

# hrs of platform (IAS) in district (D) assigned 
to mission area (M) at class (L) performed 
for air missions 

# hrs of platform (IDA) in district (D) assigned 
to mission area (M) at class (L) performed 
for surface missions 



EQUATIONS 
OBJ1 
OBJ 2 
OBJ 3 

AVAILAS (IAS) 
AIRHRREQ ( D , M , L) 
SURHRREQ ( D , M , L) 

PCTAIR(D,M) 

PCTSUR (D,M) 

DISSURF ( IDS , D) 
DISAIR(IDA, D) 

AREA (IAS) 

USEDSUR ( IDS , D , M , L) 

USEASUR ( IAS , M , L) 
USEDAIR (IDA, D,M, L) 
ACQ 
MOV 

SUPXDS (IDS , D) 
DEMXDS ( IDS , DP) 
SUPXDA ( I DA , D ) 
DEMXDA ( IDA, DP) 



total operating cost 

operating cost without acquisition costs 
operating cost without movement costs 
observe current force structure of IAS 
meet all AIR resource hour requirements 
meet all SURFACE resource hour 
requirements 

satistfy AIR TO SURFACE MAX PERCENTAGE 
satistfy SURFACE TO AIR MAX PERCENTAGE 
Produce enough district surface platforms 
Produce enough district air platforms 
Produce enough area surface platforms 
Satisfy max usage of district surface 
units 

Satisfy max usage of area surface units 

Satisfy max usage of district air units 

Limit total acquisition costs 

Limit total movement costs 

Do not exceed supply of ids platforms 

Meet demand for ids platforms 

Do not exceed supply of ida platforms 

Meet demand for ida platforms; 



* minimize 
OBJ1. . 

TCOST =E= SUM (IAS, COSTAS ( IAS )* (OXAS ( IAS ) + NXAS(IAS))) 

+ SUM ( ( IDS , D) , COSTDS ( IDS ) * (OXDS ( IDS , D) + NXDS ( IDS , D) ) ) 
+ SUM ( (IDA, D) , COSTDA ( IDA) * (OXDA (IDA, D) + NXDA ( IDA, D) ) ) 
+ SUM (IAS, PENALTY1*NXAS (IAS) ) 

+ SUM ( (IDS, D), PENALTY1*NXDS (IDS , D) ) 

+ SUM ( (IDA, D) , PENALTY1*NXDA(IDA, D) ) 

+ PENALTY 2 * SUM ( (IDA, D, DP) $ (ORD (D) NE ORD(DP)), 

MA (IDA, D, DP) ) 

+ PENALTY2 * SUM ( (IDS , D, DP) $ (ORD (D) NE ORD(DP) ) , 

MS (IDS, D, DP) ) ; 
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0BJ2$ (pick eq 2).. 

TCOST =E= SUM (IAS, COSTAS (IAS) * (OXAS (IAS) + NXAS(IAS))) 

+ SUM ( ( IDS , D) , COSTDS (IDS) * (OXDS (IDS ,D) + NXDS (IDS , D) ) ) 
+ SUM ( ( IDA, D) , COSTDA (IDA) * (OXDA (IDA, D) + NXDA (IDA, D) ) ) 
+ SUM ( IAS , PENALTY 1*NXAS (IAS) ) 

+ SUM ( (IDS , D) , PENALTY1*NXDS (IDS , D) ) 

+ SUM ( (IDA, D) , PENALTY1*NXDA(IDA,D) ) ; 

OBJ3$(pick eq 3).. 

TCOST =E= SUM (IAS, COSTAS (IAS) * (OXAS (IAS) + NXAS(IAS))) 
+ SUM ( (IDS , D) , COSTDS ( IDS ) * (OXDS ( IDS , D) + NXDS (IDS , D) ) ) 
+ SUM ( (IDA, D) , COSTDA (IDA)* (OXDA ( IDA, D) + NXDA ( IDA, D) ) ) 
+ PENALTY2 * SUM( (IDA, D, DP) $ (ORD(D) NE ORD(DP)), 

MA (IDA, D, DP) ) 

+ PENALTY2 * SUM ( ( IDS , D, DP) $ (ORD ( D) NE ORD(DP) ) , 

MS (IDS, D, DP) ) ; 



* subject to 

AVAILAS (IAS) . . OXAS (IAS) =L=SUM(D, CURRENT ( IAS , D) ) ; 

AIRHRREQ (D,M, L) $ (HRAIR(D,M, L) GT 0).. 

SUM ( IDA$ (ULDA ( IDA, M, L) GT 0), 

YDA (IDA, D,M, L) ) + 

SUM(IDS$ (ULDS (IDS ,M, L) GT 0 AND SURPCT(M) GT 0) , 
ZDS (IDS , D,M, L) ) + 

SUM ( IAS$ (ULAS ( IAS , M , L) GT 0 AND SURPCT(M) GT 0) , 
ZAS ( IAS , D , M , L) ) =G= HRAIR(D,M,L) ; 

SURHRREQ ( D , M , L) $ (HRSURF (D,M, L) GT 0).. 

SUM ( IDA$ (ULDA (IDA, M, L) GT 0 AND AIRPCT(M) GT 0) , 
ZDA (IDA, D,M, L) ) + 

SUM(IDS$ (ULDS (IDS, M,L) GT 0), 

YDS (IDS, D,M, L) ) + 

SUM ( IAS $ (ULAS ( IAS, M,L) GT 0), 

YAS ( IAS , D , M , L) ) =G= HRSURF (D,M, L) ; 

PCTAIR ( D, M) $ (AIRPCT (M) GT 0).. 

SUM ( (L, IDA) $ (ULDA (IDA, M,L) GT 0), 

ZDA(IDA,D,M,L) ) =L= 

AIRPCT (M) * SUM (L, HRSURF (D,M, L) ) ; 

PCTSUR(D,M) $ (SURPCT (M) GT 0).. 

SUM ( ( IAS , L) $ (ULAS ( IAS , M , L) GT 0), 

ZAS (IAS , D,M, L) ) + 

SUM ( (IDS , L) $ (ULDS (IDS , M , L) GT 0), 

ZDS (IDS , D, M, L) ) =L= 

SURPCT (M) * SUM ( L, HRAIR(D,M,L) ) ; 



DISSURF (IDS , D) . . 

SUM ( (M, L) $ (ULDS (IDS ,M, L) GT 0), 

YDS (IDS , D,M, L) + ZDS (IDS , D, M, L) ) 
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=L= RHDS(IDS) * (OXDS (IDS, D) +NXDS (IDS, D) ) ; 



DISAIR( IDA, D) . . 

SUM( (M,L) $ (ULDA(IDA,M,L) GT 0) , 

YDA(IDA, D,M, L) + ZDA(IDA,D,M, L) ) 

=L= RHDA(IDA) * (OXDA ( IDA, D) +NXDA (IDA, D) ) ; 

AREA ( IAS ) . . 

SUM ( (D,M, L) $ (ULAS (IAS,M, L) GT 0), 

YAS ( IAS , D,M, L) + ZAS (IAS , D,M, L) ) 

=L= RHAS(IAS) * ( OXAS ( IAS ) +NXAS ( IAS ) ) ; 



USEDSUR ( IDS , D, M, L) $ (ULDS ( IDS , M, L) GT 0).. 

YDS ( IDS , D,M, L). + ZDS (IDS , D, M, L) 

=L= (OXDS (IDS, D)+NXDS (IDS, D) ) 

* RHDS (IDS) * ULDS (IDS ,M, L) ; 



USEASUR ( IAS, M,L)$ (ULAS ( IAS, M,L) GT 0).. 

SUM (D, YAS ( IAS , D , M , L) + ZAS ( IAS , D,M, L) ) 

=L= RHAS (IAS) * ULAS ( IAS , M , L) * 

( OXAS ( I AS ) +NXAS ( IAS ) ) ; 
USEDAIR(IDA,D,M,L) $ (ULDA (IDA,M, L) GT 0).. 

YDA(IDA, D,M, L) + ZDA(IDA, D,M, L) 

=L= RHDA(IDA) * (OXDA(IDA, D) +NXDA (IDA, D) ) 
* ULDA (IDA, M, L) ; 



SUPXDS ( I DS , D ) . . 
DEMXDS (IDS, DP) . . 
SUPXDA (IDA, D) . . 
DEMXDA (IDA, DP) . . 



SUM (DP, MS ( IDS , D, DP) ) =L= CURRENT ( IDS , D) 
SUM (D, MS (IDS, D, DP) ) =E= OXDS ( IDS , DP) ; 
SUM ( DP , MA (IDA, D, DP) ) =L= CURRENT (IDA, D) 
SUM(D, MA (IDA, D, DP) ) =E= OXDA (IDA, DP) ; 



ACQ$ (ACQLIM EQ 1).. 

PENALTY 1 * (SUM ( (D, IDA) , NXDA (IDA, D) ) + 

SUM ( (D, IDS) , NXDS (IDS , D) ) + 

SUM (IAS, NXAS (IAS) ) ) =L= LIMACQ; 

MOV$ (MOVLIM EQ 1) . . 

PENALTY 2 * SUM ( ( IDA, D , DP) $ (ORD (D) NE ORD(DP)), 

MA(IDA,D,DP) ) + 

PENALTY 2 * SUM( (IDS, D, DP) $ (ORD(D) NE ORD(DP)), 

MS (IDS, D, DP) ) =L= LIMMOV; 



YAS . FX ( IAS , INDP, M, L) $ (CARD ( INDP) GT 0) = 0.0; 
YAS . FX ( IAS , D, M, L) $ (NOIAS EQ 0) = 0.0; 

YDS . FX ( IDS , D, M , L) $ (NOIDS EQ 0) = 0.0; 

YDA. FX(IDA, D,M, L) $ (NOIDA EQ 0) = 0.0; 

YAS . FX ( IAS , D,M, L) $ (ULAS (IAS , M, L) EQ 0) = 0.0; 
YDS . FX (IDS , D,M , L) $ (ULDS ( IDS, M, L) EQ 0) = 0.0; 
YDA. FX (IDA, D,M, L) $ (ULDA (IDA, M, L) EQ 0) = 0.0; 
ZDA . FX ( IDA, D,M, L) $ (ULDA ( IDA,M, L) EQ 0) = 0.0; 
ZAS. FX (IAS, D, M, L) $ (ULAS (IAS, M,L) EQ 0) = 0.0; 
ZDS . FX (IDS , D, M, L) $ (ULDS (IDS, M, L) EQ 0) = 0.0; 
ZDA. FX (IDA, D,M, L) $ (AIRPCT (M) EQ 0) = 0.0; 
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ZAS . FX ( IAS , D, M, L) $ ( SURPCT (M) EQ 0) = 0.0; 

ZDS . FX ( IDS , D, M, L) $ ( SURPCT (M) EQ 0) = 0.0; 

MODEL FLEETMIX /ALL/; 

SOLVE FLEETMIX USING RMIP MINIMIZING TCOST; 

CNT = SUM ( ( IDS , D) , OXDS . L(IDS , D) + NXDS . L ( IDS , D) ) + 
SUM ( (IDA, D) , OXDA. L(IDA, D) + NXDA. L ( IDA, D) ) + 
SUM ( IAS , OXAS.L(IAS) + NXAS . L ( IAS ) ) ; 

EQUATIONS 

TOTAL minimum total number of platforms needed; 

TOTAL.. SUM ( (IDS , D) , OXDS(IDS,D) + NXDS ( IDS , D) ) + 

SUM ( (IDA, D) , OXDA (I DA, D) + NXDA ( I DA, D)) + 

SUM (IAS, OXAS (IAS) + NXAS (IAS)) =G= FLOOR ( CNT ) ; 
MODEL FLEET /ALL/ 

SOLVE FLEET USING MIP MINIMIZING TCOST; 



* Reports 

* print the optimal objective value and solution 

OPTION DECIMALS = 5; 

DISPLAY TCOST. L, CNT; 

PARAMETER COST (*,*); 



COST ( IDS, D) 

COST (IDS, 'TYPE TOTAL') 

COST (IDS, 'UNIT COST') 

COST (IDS, 'TOTAL COST') 

COST (IDA, D) 

COST (IDA, 'TYPE TOTAL') 

COST (IDA, 'UNIT COST') 

COST (IDA, 'TOTAL COST') 

COST (IAS, 'TYPE TOTAL') 

COST (IAS, 'UNIT COST') 

COST (IAS, 'TOTAL COST') 

COST ( ' TOTAL ' , ' TOTAL COST ' ) 



OPTION DECIMALS = 2; 
DISPLAY COST; 



= OXDS . L (IDS , D) ; 

= SUM (D, OXDS . L (IDS , D) ) ; 

= COSTDS (IDS) * 1000000; 

= COSTDS (IDS) * 1000000 * 

SUM ( D, OXDS . L(IDS , D) ) ; 

= OXDA. L( IDA, D) ; 

= SUM (D, OXDA. L( IDA, D) ) ; 

= COSTDA (IDA) * 1000000; 

= COSTDA (IDA) * 1000000 * 

SUM (D, OXDA. L(IDA, D) ) ; 

= OXAS . L ( IAS ) ; 

= COSTAS (IAS) * 1000000; 

= COSTAS (IAS) *1000000* 

OXAS . L ( IAS ) ; 

= SUM (IDS, COST (IDS, 'TOTAL 
COST ' ) ) + 

SUM (IDA, COST (IDA, 'TOTAL 

COST ' ) ) + 

SUM ( IAS , COST ( IAS , • TOTAL 

COST ' ) ) ; 



PARAMETER SHORTFALL (*,* ) ACQUISITIONS RECOMMENDED; 

SHORTFALL ( I DS,D) = NXDS . L(IDS , D) ; 

SHORTFALL (IDS, 'TYPE TOTAL' ) = SUM(D, NXDS . L(IDS , D) ) ; 
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SHORT FALL (IDS, 'UNIT COST') 
SHORTFALL ( IDS , 'ACQ COST') 
SHORTFALL (IDS,' TOTAL COST ' ) 

SHORTFALL (I DA, D) 

SHORTFALL (IDA, 'TYPE TOTAL' ) 
SHORTFALL (IDA, 'UNIT COST') 
SHORTFALL (IDA, 'ACQ COST') 
SHORTFALL ( IDA , ' TOTAL COST ' ) 

SHORTFALL (IAS, 'TYPE TOTAL') 
SHORTFALL ( IAS , ' UNIT COST ' ) 
SHORTFALL ( IAS , ' ACQ COST ' ) 
SHORTFALL ( IAS , ' TOTAL COST ' ) 

SHORTFALL ( 'TOTAL' , 'TOTAL COST' ) = 



COSTDS (IDS) *1000000; 
PENALTY1*1000000; 

(COSTDS (IDS) +PENALTY1) *1000000 

* SUM(D, NXDS . L (IDS , D) ) ; 
NXDA.L(IDA,D) ; 

SUM(D, NXDA. L(IDA, D) ) ; 

COSTDA (IDA) *1000000; 
PENALTY1* 10000 00 ; 

( COSTDA ( IDA) +PENALTY1) *1000000 

* SUM(D, NXDA. L (IDA, D) ) ; 
NXAS.L(IAS) ; 

COSTAS (IAS) *1000000; 
PENALTY1*1000000 ; 

(COSTAS (IAS) +PENALTY1) *1000000 

* NXAS.L(IAS); 

SUM ( IDS , SHORTFALL ( IDS , ' TOTAL 
COST ' ) ) 

+ SUM (IDA, SHORTFALL (IDA,' TOTAL 
COST ' ) ) 

+ SUM ( IAS , SHORTFALL ( IAS , ' TOTAL 
COST ' ) ) ; 



OPTION DECIMALS = 0; 

DISPLAY $(SHORTFALL( 'TOTAL' , 'TOTAL COST') GT 0) SHORTFALL; 

PARAMETER MOVEMENT (*, D, DP) PLATFORM MOVEMENT FROM (ROW) TO 

(COLUMN) ; 

MOVEMENT (IDS, D, DP) = MS . L(IDS , D, DP) ; 

MOVEMENT ( IDA, D, DP) = MA. L(IDA, D, DP) ; 

OPTION MOVEMENT : 2 ; 1 : 1 ; 

DISPLAY MOVEMENT; 



PARAMETER ALLOCATE (D, P,M) , TALL(*,M) ; 

ALLOCATE (D, I DS,M) = SUM(L, ZDS . L (IDS , D,M, L) +YDS . L(IDS , D, M, L) ) 

* 100 . 0 ; 

ALLOCATE (D, IDA, M) = SUM(L, ZDA. L(IDA, D,M, L) +YDA. L(IDA, D,M, L) ) 

* 100 . 0 ; 

ALLOCATE (D, I AS, M) = SUM (L, ZAS . L (IAS , D,M, L) +YAS . L ( IAS , D, M, L) ) 

* 100 . 0 ; 

OPTION ALLOCATE : 2 : 1 : 1 ; DISPLAY ALLOCATE; 

TALL(P,M) = SUM(D, ALLOCATE (D, P,M) ); 

TALL ( 'Total' ,M) =SUM( (D,P), ALLOCATE (D, P,M) ); 

DISPLAY TALL; 
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APPENDIX C 



GAMS OUTPUT 



526 VARIABLE TCOST.L - 172.11000 TOTAL OPERATING COST OF FLEET PER YEAR 

PARAMETER CNT - 147.91996 PLATFORM COUNT 



547 PARAMETER COST 





FIRST 


FIFTH 


SEVENTH 


EIGHTH 


HU25 


3.00 




4.00 




HH65 


5.00 


3.00 


16.00 


8.00 


HH3 


4.00 


2.00 


8.00 




C130 




2.00 


3.00 




WPB82 


8.00 


6.00 


8.00 


5.00 


WPB110 


2.00 


2.00 


16.00 


2.00 


WTGB 


3.00 








WLB180 




1.00 


2.00 


2.00 


WLM157 


1.00 


2.00 






WLM133 


2.00 






1.00 



WHEC 

WMEC270 

WMEC210 

TOTAL 

583 PARAMETER MOVEMENT 



NINTH 


TYPE TOTAL 


UNIT COST 


TOTAL COST 




7.00 


1640000.00 


11480000.00 


2.00 


34.00 


1110000.00 


37740000.00 


1.00 


15.00 


1470000.00 


22050000.00 




5.00 


2220000.00 


11100000.00 




27.00 


370000.00 


9990000.00 




22.00 


590000.00 


12980000.00 


4.00 


7.00 


650000.00 


4550000.00 


2.00 


7.00 


1580000.00 


11060000.00 




3.00 


1130000.00 


3390000.00 


2.00 


5.00 


990000.00 

4280000.00 


4950000.00 




8.00 


2690000.00 


21520000.00 




10.00 


2130000.00 


21300000.00 
1 . 721100E+8 



PLATFORM MOVEMENT FROM (ROW) TO (COLUMN) 



INDEX 1 = HU25 



FIRST 

FIRST 3.00 

SEVENTH 

INDEX 1 - HH65 



SEVENTH 

4.00 



FIRST FIFTH SEVENTH EIGHTH NINTH 



FIRST 

FIFTH 

SEVENTH 

EIGHTH 

NINTH 



5.00 



3.00 



16.00 



8.00 



2.00 
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583 PARAMETER MOVEMENT 



PLATFORM MOVEMENT FROM (ROW) TO (COLUMN) 



INDEX 1 ■= HH3 



FIRST 

FIRST 4.00 

FIFTH 

SEVENTH 

NINTH 



INDEX 1 - C130 

FIFTH 

FIFTH 2.00 

SEVENTH 

INDEX 1 - WPB82 

FIRST 

FIRST 8.00 

FIFTH 

SEVENTH 

EIGHTH 

INDEX 1 - WPB110 

FIRST 

FIRST 2.00 

FIFTH 

SEVENTH 

EIGHTH 

INDEX 1 - WTGB 

FIRST 

FIRST 3.00 

NINTH 

INDEX 1 - WLB180 

FIFTH 

FIFTH 1.00 

SEVENTH 

EIGHTH 

NINTH 



FIFTH SEVENTH 

2.00 

8.00 

SEVENTH 

3.00 

FIFTH SEVENTH 

6.00 

8.00 

FIFTH SEVENTH 

2.00 

16.00 

NINTH 

4.00 

SEVENTH EIGHTH 

2.00 

2.00 



NINTH 



1.00 



EIGHTH 



5.00 



EIGHTH 



2.00 



NINTH 



2.00 
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583 PARAMETER MOVEMENT 



PLATFORM MOVEMENT FROM (ROW) TO (COLUMN) 



INDEX 1 


- WLM157 






FIRST 


FIFTH 


FIRST 

FIFTH 


1.00 


2.00 


INDEX 1 


- WLM133 






FIRST 


SEVENTH 


FIRST 

SEVENTH 

EIGHTH 


2.00 


1.00 



EIGHTH 



2.00 



591 PARAMETER ALLOCATE 



INDEX 1 = FIRST 





SAR 


ELT 


ATN 


HU25 


939.50 


287.50 


200.00 


HH65 


768.50 


2412.50 


44.00 


HH3 


160.00 


1336.00 




WPB82 


3330.00 


9320.00 




WPB110 


705.00 


2400.00 




WTGB 


2030.00 


2070.00 




WLM157 




40.00 


3560.00 


WLM133 




336.00 


5565.00 


WMEC270 




850.00 





INDEX 1 


= FIFTH 








SAR 


ELT 


ATN 


HH65 




1275.00 


171.00 


HH3 

C130 


260.00 

1600.00 


995.00 




WPB82 

WPB110 

WLB180 


1580.00 


7990.00 

2400.00 


3876.50 


WLM157 




562.50 


7437.50 


WMEC270 

WMEC210 


681.00 


1892.50 





IOP 


MILOP 


MSAFE 


MEP 


4.00 






528.00 


953.00 


1085.00 


1133.00 

895.00 


15.00 

550.00 


481.00 


770.00 




548.00 


IOP 


MILOP 


MSAFE 


V 

MEP 


1.00 


221.00 

70.00 


18.00 


249.00 




1600.00 

123.50 


59.00 


271.00 




700.00 

641.50 
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591 PARAMETER ALLOCATE 



INDEX 1 


= SEVENTH 














SAR 


ELT 


ATN 


MILOP 


MSAFE 


MEP 


HU25 


2221.00 


900.00 


79.00 








HH65 




10320.00 










HH3 


2296.00 


2800.00 






16.00 


170.00 


C130 


123.00 


2048.00 




229.00 






WPB82 


5200.00 


7848.00 






66.00 


86.00 


WPB110 


255.00 


30262.00 




1483.00 






WLB180 




669.00 


7331.00 








WLM133 




231.00 


3619.00 








WMEC270 




27171.50 










WMEC210 




38868.50 










INDEX 1 


- EIGHTH 














SAR 


ELT 


ATN 


MILOP 


MSAFE 


MEP 


HH65 


1670.00 


3020.00 


63.00 


161.00 


14.00 


206.00 


WPB82 


2420.00 


5830.00 










WPB110 




3438.00 




512.00 




50.00 


WLB180 




1058.00 


6579.00 


300.00 


63.00 




WLM133 




509.00 


7191.00 








WMEC270 




3265.00 










WMEC210 




3690.00 










INDEX 1 


- NINTH 














SAR 


ELT 


ATN IOP 


MILOP 


MSAFE 


MEP 


HH65 


354.00 


42.00 169 


.00 73.00 


36.00 


15.00 


55.00 


HH3 


700.00 












WTGB 


140.00 


163.00 


2030.00 


35.00 


136.00 




WLB180 




126.00 6069 


.00 1800.00 






5.00 



594 


PARAMETER 


TALL 














SAR 


ELT 


ATN 


IOP 


MILOP 


MSAFE 


MEP 


HU25 


3160 


1187 


279 


4 






528 


HH65 


2792 


17069 


447 


74 


418 


47 


510 


HH3 


3416 


5131 






1203 


31 


170 


C130 


1723 


2048 






229 






WPB82 


12530 


30988 








675 


357 


WPB110 


960 


38500 






4490 




50 


WTGB 


2170 


2233 




2983 


1120 


136 




WLB180 




1853 


23855 


1800 


423 


63 


5 


WLM157 




1002 


10997 










WLM133 




1076 


16375 


481 


770 




548 


WMEC270 


681 


33179 






700 






WMEC210 




42558 






641 






TOTAL 


27433 


176826 


51954 


5342 


9995 


952 


2168 
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APPENDIX D. FORTRAN LISTING 



PROGRAM USCG 

******************************************************************* 



* * 

* PROGRAM: USER INTERFACE FOR COAST GUARD FORCE STRUCTURE * 

* OPTIMIZATION MODEL * 

* * 

* AUTHOR : LT JOHN E TOMKO, USCG * 

* * 

* WRITTEN: 13 MAY 1991 * 

* * 

* LAST MODIFIED: 01 SEPTEMBER 1991 * 

* * 

* LANGUAGE : FORTRAN * 

* * 

* PROBLEM STATEMENT: A MENU-DRIVEN FRONT END USER INTERFACE FOR * 

* MANIPULATING DATASETS CONTAINED IN SEPERATE FILES TO PRODUCE * 

* A GAMS FORMAT FOR DATA TO BE USED IN AN OPTIMIZATION MODEL. * 

* * 



******************************************************************* 



* EXTERNAL FILES: 

* INPUT : 



* 

* 

* 

* 

* 

* 

* 

* 

* 

* 



AMISS DATA A1 
APLAT DATA Al 
ADI ST DATA Al 
AHOURS DATA Al 
ATABLE DATA Al 
PMISS DATA Al 
PPLAT DATA Al 
PD I ST DATA Al 
PHOURS DATA Al 
PTABLE DATA Al 



* OUTPUT : 



* 

* 

ATLANTIC AREA MISSION DATA * 

LANT AREA PLATFORM NAMES & DATA * 

ATLANTIC AREA DISTRICT DATA * 

LANT AREA RESOURCE HOUR REQUIREMENTS * 

LANT AREA PARAMETER DATA * 

PACIFIC AREA MISSION DATA * 

PAC AREA PLATFORM NAMES & DATA * 

PACIFIC AREA DISTRICT DATA * 

PAC AREA RESOURCE HOUR REQUIREMENTS * 

PAC AREA PARAMETER DATA * 

* 



* GAMS DATA Al : LISTING OF GAMS DATA FOR $ INCLUDE * 

******************************************************************* 



INTEGER I, J, K, L, NUMDIST, NUMPLAT(A) , NUMISS , HRS(24) 
INTEGER HOURS (7, 10, 6) , QTY(24,7), CATB(24,3), CATA(24,10) 
INTEGER INDP , LOCINDP(7) 

REAL VALA(4) , VALB(4) , ATOS(IO), STOA(IO) , COST(24) , PCNT 
REAL MODEL(7) 

CHARACTER MISSIONS (10)*5 , DISTRICT(7)*10 , PLATFORM (24) *8 
CHARACTER WORDA(4)*18, WORDB(4)*18, TYPE(24)*3, AREA* 8 
LOGICAL GAMS 
GAMS = .FALSE. 

CALL INITIAL(NUMDIST , NUMPLAT , NUMISS , HRS , HOURS , COST , CATA , 
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+ CATB , VALA , VALB , ATOS , STOA , QTY) 

CALL GETDATA (NUMDI ST , NUMPLAT , NUMI S S , HRS , HOURS .COST, MISSIONS, 

+ DISTRICT , PLATFORM , TYPE , VALA , CATA , CATB , VALB , ATOS , 

+ STOA , AREA , QTY , WORDA , WORDB , INDP , LOCINDP ) 

CALL MAINMENU(NUMDIST , NUMPLAT , NUMI SS , HRS , HOURS , COST , MISSIONS , 

+ DISTRICT, PLATFORM , TYPE , CATA , CATB , VALA , VALB , ATOS , 

+ STOA , AREA , MODEL , QTY , WORDA , WORDB , GAMS , INDP , 

+ LOCINDP) 

CALL SAVEDATA(NUMDIST , NUMPLAT , NUMISS , HRS , HOURS , COST , MISSIONS , 

+ DISTRICT , PLATFORM , TYPE , VALA , VALB , ATOS , STOA , AREA , 

+ QTY , WORDA , WORDB , CATA , CATB , INDP , LOCINDP) 

IF (GAMS) THEN 
WRITE(* , 11) 

ELSE 

WRITE(* , 22) 

ENDIF 

11 FORMAT(/ , ' TO RUN MODEL TYPE: GAMS FLEET') 

22 FORMAT(/ , ' NO GAMS MODEL WAS GENERATED DURING THIS SESSION') 
STOP 
END 



SUBROUTINE GETDATA (NUMDI ST , NUMPLAT , NUMISS , HRS , HOURS , COST , 

+ MI SS IONS , DI STRICT , PLATFORM , TYPE , VALA , CATA , 

+ CATB , VALB , ATOS , STOA , AREA , QTY , WORDA , WORDB , 

+ INDP, LOCINDP) 

********************************************************************** 



* * 

* SUBROUTINE: GET DATA FROM EXTERNAL FILES * 

* * 

* AUTHOR : LT J.E. TOMKO , USCG * 

* * 

* WRITTEN : 19 MAY 1991 * 

* * 

* MODIFIED : 09 JUNE 1991 * 

* * 
********************************************************************** 

* PURPOSE: READ THE EXTERNAL FILES INTO ARRAYS FOR USE WITH THE * 

* INTERFACE. * 

* * 



********************************************************************** 
INTEGER I, J, K, L, NUMDIST, NUMPLAT(A) , NUMISS, HRS(24) 

INTEGER HOURS (7 ,10,6) , QTY(24,7), CATB(24,3), CATA(24,10) 

INTEGER INDP, LOCINDP(7) 

REAL VALA(4) , VALB(4) , ATOS(IO), STOA(IO) , COST(24) 

CHARACTER MISSIONS ( 10)*5 , DISTRICT(7)*10 , PLATFORM (24) *8 
CHARACTER WORDA(4)*18, WORDB(4)*18, TYPE(24)*3, T*l, AREA* 8 
LOGICAL WRONG 
WRONG - .TRUE. 

DO 5 WHILE (WRONG) 

WRITE(* , 111) 

READ(* , 99) T 
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.FALSE. 



IF ((T .EQ. 'P') .OR. (T .EQ. 'A')) WRONG - 
5 ENDDO 

OPEN(21,FILE-'/'//T//'MISS DATA Al') 

OPEN (22, FILE- '/'//I //' PLAT DATA Al') 
OPEN(23,FILE-'/'//T//'DIST DATA Al') 

OPEN(24 , FILE—'/' //T//'HOURS DATA Al') 

OPEN(26 , FILE-'/' //T//' TABLE DATA Al') 

PRINT* ,' RETRIEVING MISSION DATA...' 

READ(21 , 11) NUMISS 
DO 10 I - 1, NUMISS 

READ(21 , 22) MISSIONS(I) 

10 CONTINUE 

PRINT*, 'RETRIEVING PLATFORM DATA. . . ' 

DO 15 J - 2,4 
NUMPLAT(J) - 0 
15 CONTINUE 

READ(22 , 11) NUMPLAT ( 1 ) 

DO 20 I = 1, NUMPLAT(l) 

READ(22 ,44) PLATFORM(I) ,TYPE(I) ,HRS(I) ,COST(I) 
IF (TYPE (I) .EQ. 'IDS') THEN 
NUMPLAT ( 2 ) = NUMPLAT(2) + 1 
ELSEIF (TYPE(I) .EQ. 'IDA') THEN 
NUMPLAT ( 3 ) - NUMPLAT ( 3 ) + 1 
ELSE 

NUMPLAT (4) - NUMPLAT (4) + 1 
END IF 

20 CONTINUE 

PRINT*, 'RETRIEVING DISTRICT DATA...' 

READ(23 , 122) AREA 
READ(23 , 11) NUMDIST 
DO 30 I - 1, NUMDIST 

READ(23 , 33) DISTRICT(I) 

30 CONTINUE 

READ(23 , 11) INDP 
DO 35 I - 1, INDP 
READ(23 , 11) LOCINDP(I) 

35 CONTINUE 

PRINT*, 'RETRIEVING PARAMETER DATA...' 

DO 40 I - 1, NUMDIST 
DO 50 J - 1, NUMISS 

READ (24 ,55) (HOURS (I , J ,K) ,K-1 , 6) 

50 CONTINUE 

40 CONTINUE 

READ (26 ,77) (VALA(I), 1-1,4) 

DO 80 I - 1,4 

READ(26 , 133) WORDA(I) 

80 CONTINUE 

DO 90 I - 1, NUMPLAT(l) 

READ(26 , 88) (CATB(I,J), J-1,3) 

90 CONTINUE 

READ (26 ,77) (VALB(I) , 1-1,4) 

DO 100 1-1,4 
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READ(26 , 133) WORDB(I) 

100 CONTINUE 

DO 110 I - 1, NUMPLAT ( 1 ) 

READ(26 ,88) (CATA(I,J), J-l, NUMISS) 

110 CONTINUE 

DO 130 I - 1, NUMPLAT (1) 

READ(26 , 88) (QTY(I,J), J-l.NUMDIST) 

130 CONTINUE 

READ(26 ,77) (STOA(I), 1-1, NUMISS) 

READ(26 ,77) (ATOS(I), 1-1, NUMISS) 

11 FORMAT (13) 

22 FORMAT ( IX , A5 ) 

33 FORMAT ( IX, A10) 

44 FORMAT ( IX , A8 , A3 , IX , 14 , IX , F7 . 4 ) 

55 FORMAT ( IX , 6 ( I 6 , IX) ) 

77 FORMAT ( 10 ( F4 . 2 , IX) ) 

88 FORMAT (7(14, IX)) 

99 FORMAT (Al) 

111 FORMAT ( ' PLEASE CHOOSE WHICH AREA DATA YOU WISH TO USE:',/, 

+ ' (A)TLANTIC' ,/, ' (P)ACIFIC' ,/, ' ENTER A OR P:') 

122 FORMAT ( IX, A8) 

133 FORMAT ( IX, Al 8) 

RETURN 

END 



SUBROUTINE SAVEDATA(NUMDIST , NUMPLAT , NUMISS , HRS , HOURS , COST , 

+ MISSIONS, DISTRICT, PLATFORM, TYPE, VALA.VALB, 

+ ATOS , STOA , AREA , QTY , WORDA , WORDB , CATA , CATB , 

+ INDP , LOCINDP) 

********************************************************************** 



* * 

* SUBROUTINE: SAVE CHANGES TO DATABASES * 

* * 

* AUTHOR : LT J.E. TOMKO, USCG * 

* * 

* WRITTEN : 19 MAY 1991 * 

* * 

* MODIFIED : 09 JUNE 1991 * 

* * 
********************************************************************** 

* PURPOSE: IF SO DESIRED, STORE CHANGES TO THE DATABASE IN THE * 

* ORIGINAL INPUT FILES FOR FUTURE USE. * 

* * 



********************************************************************** 
INTEGER I, J, K, L, NUMDIST, NUMPLAT (4) , NUMISS, HRS(24) 

INTEGER HOURS (7, 10, 6) , QTY(24,7), CATB(24,3), CATA(24,10) 

INTEGER INDP, LOCINDP(7) 

REAL VALA(4) , VALB(4) , ATOS(IO), STOA(IO) , COST(24) 

CHARACTER MISSIONS(10)*5, DISTRICT(7)*10, PLATFORM (24) *8 
CHARACTER W0RDA(4)*18, WORDB(4)*18, TYPE(24)*3, ANS*1, AREA* 8 
REWIND (21) 
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10 

20 

30 

35 

50 

40 

70 

80 

90 

100 

110 



REWIND(22) 

REWIND (23) 

REWIND (24) 

REWIND (26) 

WRITE(* , 11) 

READ(*,22) ANS 
IF (ANS .EQ. 'S') THEN 

PRINT*, 'SAVING MISSION DATA...' 

WRITE(21 , 33) NUMISS 
DO 10 I - 1, NUMISS 

WRITE(21 , 77) MISSIONS(I) 

CONTINUE 

PRINT* ,' SAVING PLATFORM DATA...' 

WRITE(22 , 33) NUMPLAT ( 1 ) 

DO 20 I - 1, NUMPLAT ( 1 ) 

WRITE(22 , 99) PLATFORM(I) ,TYPE(I) ,HRS(I) ,COST(I) 
CONTINUE 

PRINT* ,' SAVING DISTRICT DATA...' 

WRITE(23 , 133) AREA 
WRITE(23 , 33) NUMDIST 
DO 30 I - 1, NUMDIST 

WRITE(23 , 88) DISTRICT(I) 

CONTINUE 

WRITE(23 , 33) INDP 
DO 35 I - 1, INDP 

WRITE(23 , 33) LOCINDP(I) 

CONTINUE 

PRINT* , ' SAVING RESOURCE HOUR DATA...' 

DO 40 I = 1, NUMDIST 
DO 50 J - 1, NUMISS 

WRITE(24, 55) (HOURS(I , J ,K) ,K-1 , 6) 

CONTINUE 

CONTINUE 

PRINT*, 'SAVING PARAMETER DATA. . . ' 

WRITE(26 , 111) (VALA(I) , 1=1,4) 

DO 70 I = 1,4 

WRITE(26 , 144) WORDA(I) 

CONTINUE 

DO 80 I = 1, NUMPLAT(l) 

WRITE(26 , 122) (CATB(I,J), J=l,3) 

CONTINUE 

WRITE(26 , 111) (VALB(I) , 1=1,4) 

DO 90 I = 1,4 

WRITE(26 , 144) WORDB(I) 

CONTINUE 

DO 100 I = 1, NUMPLAT ( 1 ) 

WRITE(26 , 122) (CATA(I,J), J=l, NUMISS) 

CONTINUE 

DO 110 I = 1, NUMPLAT(l) 

WRITE(26 , 122) (QTY(I.J), J-l, NUMDIST) 

CONTINUE 

WRITE(26 , 111) (STOA(I) , 1=1, NUMISS) 



73 



WRITE(26 , 111) (ATOS(I) , 1=1, NUMISS) 

ELSE 

WRITE(*,44) 

ENDIF 

CLOSE(21) 

CLOSE(22) 

CLOSE(23) 

CLOSE(24) 

CLOSE(26) 

II FORMAT ( ' ’/EXITING PROGRAM.',/.' YOU MAY EITHER’ (S)AVE 

+ 'CHANGES MADE DURING THIS SESSION' , 5X/OR' ,'// (E)XIT ', 

+ 'WITHOUT CHANGING THE MASTER DATABASE.’,/, 

+ ' PLEASE ENTER S OR E:') 

22 FORMAT (Al) 

33 FORMAT (13) 

44 FORMAT ( ' ’/MASTER DATABASE NOT MODIFIED BY CHANGES MADE DURING ', 
+ 'THIS SESSION.') 

55 FORMAT ( IX , 6 ( I 6 , IX) ) 

77 FORMAT ( IX , A5 ) 

88 FORMAT ( IX, AlO) 

99 F0RMAT(1X, A8 , A3 , IX, 14 , IX, F7 . 4) 

III FORMAT ( 10 (F4. 2, IX)) 

122 FORMAT (7(14, IX) ) 

133 FORMAT ( IX , A8 ) 

144 FORMAT ( ’ ' ,A18) 

RETURN 

END 



SUBROUTINE MAINMENU(NUMDIST , NUMPLAT , NUMISS , HRS , HOURS , COST , 

+ MISSIONS, DISTRICT, PLATFORM, TYPE, CATA, 

+ CATB , VALA , VALB , ATOS , STOA , AREA , MODEL , QTY , 

+ WORDA , WORDB , GAMS , PCNT , INDP , LOCINDP) 

********************************************************************** 



* * 

* SUBROUTINE: MAIN MENU CONTROLLER * 

* * 

* AUTHOR : LT J.E. TOMKO , USCG * 

* * 

★ WRITTEN : 20 MAY 1991 * 

* * 

* MODIFIED : 04 JUNE 1991 * 

★ * 
********************************************************************** 

* PURPOSE: PROVIDE USER WITH MAIN MENU WHICH CAN ACCESS ALL OTHER * 

* MENUS AVAILABLE IN THE INTERFACE. * 

* * 



********************************************************************** 
INTEGER I, J, K, L, NUMDIST, NUMPLAT (4) , NUMISS, HRS(24) 

INTEGER HOURS (7 ,10,6) , QTY(24,7), CATB(24,3), CATA(24,10) 

INTEGER INDP, LOCINDP(7) , PICK, CNT 

REAL VALA(4) , VALB(4) , ATOS(IO) , STOA(IO) , COST(24) 
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REAL M0DEL(7) 

CHARACTER MISSIONS(10)*5 , DISTRICT(7)*10 , PLATFORM(24)*8 
CHARACTER WORDA(4)*18, WORDB(4)*18, TYPE(24)*3, ANS*1 , 0PT*1 
CHARACTER AREA* 8 
LOGICAL GAMS 
ANS - 'A' 

PICK - 0 
CNT - 0 

DO 10 WHILE (ANS .NE. 'E') 

WRITE(* , 11) 

READ(*,22) ANS 
IF (ANS .EQ. 'D') THEN 
PICK - 1 

CALL DEFINE (MODEL) 

ELSEIF (ANS .EQ. 'P') THEN 

CALL OUT (NUMDIST , NUMPLAT , NUMISS , HRS , HOURS , COST , CNT .MISSIONS , 

+ DI STRICT , PLATFORM , TYPE , CATA , CATB , VALA , VALB , ATOS , 

+ STOA , QTY , PCNT , INDP , LOCINDP , WORDA , WORDB) 

ELSEIF (ANS .EQ. 'M') THEN 
OPT - 'A' 

DO 20 WHILE (OPT .NE. 'E') 

WRITE(*,33) 

READ(* , 22 ) OPT 
IF (OPT .EQ. 'P') THEN 

CALL EDITPLAT (NUMDIST , NUMPLAT , NUMISS , HRS , HOURS , COST , CATA , 
+ CATB .MISSIONS, DISTRICT , PLATFORM , TYPE , QTY , 

+ WORDA, WORDB, VALA, VALB) 

ELSEIF (OPT .EQ. 'D') THEN 

CALL EDI TDI ST (NUMDIST , NUMPLAT , NUMISS , HRS , HOURS , COST , QTY , 

+ MISSIONS , DISTRICT , PLATFORM , TYPE , AREA , INDP , 

+ LOCINDP) 

ELSEIF (OPT .EQ. 'M') THEN 

CALL EDITMISS (NUMDIST , NUMPLAT , NUMISS , HOURS , MISSIONS , 

+ DISTRICT , PLATFORM , ATOS , STOA , CATA , VALA , WORDA) 

ELSEIF (OPT .NE. 'E') THEN 
WRITE(*. 55) 

ENDIF 

20 ENDDO 

ELSEIF (ANS .EQ. 'G') THEN 
IF (PICK .NE. 0) THEN 
GAMS - .TRUE. 

CALL DOIT (NUMDIST , NUMPLAT , NUMISS , HRS , HOURS , COST , MISSIONS , 

+ DISTRICT, PLATFORM , TYPE , CATA , CATB , VALA , VALB , 

+ ATOS , STOA , QTY , PCNT , INDP , LOCINDP , MODEL) 

ELSE 

WRITE(* , 44) 

ENDIF 

ELSEIF (ANS .EQ. 'C') THEN 

CALL CHANGE (WORDA , WORDB , VALA , VALB , ATOS , STOA , NUMPLAT , 

+ NUMISS, MISSIONS, PLATFORM, PCNT) 

ELSEIF (ANS .NE. 'E') THEN 
WRITE(*,55) 
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ENDIF 

10 ENDDO 

11 FORMAT ( ' ' , 20X, 'MAIN MENU',/,' THE FOLLOWING OPTIONS ARE ', 

+ 'AVAILABLE:',/,' (D)EFINE PROBLEM' ,/, ' (P)RINT DATABASE' 

+ ,/,' (M)ODIFY DATABASE' ,/, ' (G)ENERATE GAMS FILE' ,/, 

+ ' (C)HANGE PARAMETERS' ,/, ' (E)XIT PROGRAM' ,/, 

+ ' ENTER YOUR CHOICE (D,P,M,G,C,E) : ' ) 

22 FORMAT (Al) 

33 FORMAT ( ' ',5X, 'MODIFY DATABASE MENU',/,' SELECT DATABASE TO BE ', 

+ 'MODIFIED:',/,' (P)LATFORM' ,/, ' (D)ISTRICT' ,/, ' (M)ISSION' , 

+ /,' (E)XIT MODIFY DATABASE ROUTINE' ,/, 

+ ' ENTER YOUR CHOICE (P,D,M,E):') 

44 FORMAT ( ' ','*** ERROR ***',/,' (D)EFINE PROBLEM FIRST',/) 

55 FORMAT ( ' ','*** ERROR ***',/,' INVALID OPTION',/) 

RETURN 

END 



SUBROUTINE INITI AL(NUMDIST , NUMPLAT , NUMISS , HRS , HOURS , COST , CATA , 

+ CATB , VALA , VALB , ATOS , STOA.QTY) 

********************************************************************** 



* * 

* SUBROUTINE: INITIALIZE ARRAYS * 

* * 

* AUTHOR : LT J.E. TOMKO , USCG * 

* * 

* WRITTEN : 25 MAY 1991 * 

* * 

* MODIFIED : 15 JUNE 1991 * 

* * 
********************************************************************** 

* PURPOSE: INITIALIZE VARIABLES TO PROTECT AGAINST BAD DATA BEING * 

* MISTAKENLY ENTERED IN THE ARRAYS. * 

* * 



********************************************************************** 
INTEGER I, J, K, L, NUMDIST, NUMPLAT (4 ) , NUMISS, HRS(24) 

INTEGER HOURS (7 ,10,6) , QTY(24,7), CATB(24,3), CATA(24,10) 

REAL VALA (4) , VALB(4) , ATOS(IO), STOA(IO) , COST(24) 

NUMISS - 0 
DO 10 I - 1,24 
HRS (I) - 0 
COST(I) - 0.0 
DO 20 J - 1,10 
CATA ( I , J ) - 1 
20 CONTINUE 

DO 30 J - 1,3 
CATB ( I , J ) - 1 
30 CONTINUE 

DO 40 J - 1,7 
QTY ( I , J ) - 0 
40 CONTINUE 

10 CONTINUE 



76 



DO 50 I - 1,4 
NUMPLAT ( I ) - 0 
VALA(I) - 0.0 
VALB(I) - 0.0 
50 CONTINUE 

DO 60 I - 1,10 
ATOS(I) - 0.0 
STOA(l) - 0.0 
60 CONTINUE 

NUMDIST - 0 
DO 80 I - 1,7 
DO 90 J - 1,10 
DO 100 K - 1,6 
HOURS (I , J ,K) - 0 
100 CONTINUE 

90 CONTINUE 

80 CONTINUE 

RETURN 
END 



SUBROUTINE DEFINE (MODEL) 



* * 

★ SUBROUTINE: DEFINE FLEET MIX PROBLEM * 

* * 

* AUTHOR : LT J.E. TOMKO, USCG * 

* * 

* WRITTEN : 04 JUNE 1991 * 

* * 

* MODIFIED : 01 SEPTEMBER 1991 * 

★ * 
***********************************1c***1t**it'k'k'k'k'k*'k'k'k'k'k**'k'k'k'k'k'k'k'k*'k'k*'k'k 

* PURPOSE: ALLOW THE SUER TO SPECIFY CERTAIN ASPECTS OF THE PROBLEM * 

* SUCH AS COSTS TO BE CONSIDERED. * 

* * 



********************************************************************** 
CHARACTER ANS*1 , OPT*l 
REAL MODEL(7 ) , TEMP 
LOGICAL ERROR 
ERROR - .TRUE. 

DO 10 WHILE (ERROR) 

WRITE(* ,11) 

READ (*,22) OPT 
ERROR - .FALSE. 

IF (OPT .EQ. 'A') THEN 
MODEL(l) - 1.0 
ELSEIF (OPT .EQ. 'B') THEN 
MODEL(l) - 3.0 
MODEL(2) - 0.0 
MODEL(4) - 0.0 
MODEL(6) - 0.0 
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ELSEIF (OPT .EQ. 'C') THEN 
MODEL(l) - 2.0 
MODEL(3) - 0.0 
MODEL(5) - 0.0 
MODEL(7) - 0.0 
ELSE 

WRITE(* ,33) 

ERROR - .TRUE. 

END IF 

10 ENDDO 

IF ((OPT .EQ. 'A') .OR. (OPT .EQ. 'C')) THEN 
WRITE(* , 55) 

READ(* , *) MODEL(6) 

MODEL(6) - MODEL(6) / 1000.0 
IF (MODEL(6) .LT. 0.0) MODEL(6) - 0.0 
WRITE(* , 44) 

READ(* , *) TEMP 
IF (TEMP .GE. 0.0) THEN 
MODEL(2) - 1.0 
MODEL(4) - TEMP / 1000.0 
ELSE 

MODEL(2) - 0.0 
MODEL(4) - 0.0 
END IF 
END IF 

IF ((OPT .EQ. 'A') .OR. (OPT .EQ. 'B')) THEN 
WRITE(* , 77) 

READ(* , *) MODEL(7) 

MODEL(7) - MODEL(7) / 1000.0 
IF (MODEL(7) .LT. 0.0) MODEL(7) - 0.0 
WRITE(* , 66) 

READ(* , *) TEMP 
IF (TEMP .GE. 0.0) THEN 
MODEL(3) -1.0 
MODEL(5) - TEMP / 1000.0 
ELSE 

MODEL(3) - 0.0 
MODEL(5) - 0.0 
END IF 
END IF 

WRITE(* , 88) 

IF (OPT .EQ. 'A') THEN 
WRITE(*,99) 

ELSEIF(OPT .EQ. 'B') THEN 
WRITE(* , 122) 

ELSE 

WRITE(* , 111) 

END IF 

IF (MODEL(2) .EQ. 0.0) THEN 
WRITE(* ,133) 

WRITE(* , 177) INT(MODEL(6)*1000000 . 0) 

ELSE 
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11 



22 

33 

44 



55 

66 

77 

88 

99 

111 

122 

133 

144 

155 

166 

177 

188 



WRITE(* , 144) INT(MODEL(6)*1000000.0) , INT(MODEL(4)*1000000 . 0) 
END IF 

IF (MODEL(3) .EQ. 0.0) THEN 
WRITE(* ,155) 

WRITE(* , 188) INT(MODEL( 7) *1000000 .0) 

ELSE 

WRITE(* , 166) INT(MODEL(7)*1000000 .0) , INT(M0DEL( 5) *1000000 . 0) 
ENDIF 

FORMAT (' ' , 'THE FOLLOWING MODELS ARE AVAILABLE: ' ,/, 

' (A) FULL MODEL WITH ACQUISITION AND MOVEMENT COSTS' ,/, 
' (B) MODEL WITHOUT ACQUISITION COST',/.' (C) MODEL ', 
'WITHOUT MOVEMENT COST',/,' ENTER CHOICE (A,B,C):') 
FORMAT (Al) 

FORMAT ( ' ','*** ERROR ***',/,' INVALID OPTION',/) 

,' ENTER UPPER LIMIT FOR TOTAL ACQUISITION COST:',/, 

' ENTER A (-1) FOR NO LIMIT. ENTER VALUE IN ', 

' THOUSANDS OF DOLLARS : ' ) 

,' ENTER ACQUISITION PENALTY TO USE IN THIS MODEL',/, 

' (IN THOUSANDS OF DOLLARS (1 MILLION = 1000.0):') 

,' ENTER UPPER LIMIT FOR TOTAL MOVEMENT COST:',/, 

' ENTER A (-1) FOR NO LIMIT. ENTER VALUE IN ', 

' THOUSANDS OF DOLLARS : ' ) 

,' ENTER MOVEMENT PENALTY TO USE IN THIS MODEL',/, 

' (IN THOUSANDS OF DOLLARS (1 MILLION - 1000.0):') 

, ' THE FOLLOWING MODEL HAS BEEN DEFINED:',/) 

, ' FULL MODEL WITH ACQUISITION AND MOVEMENT COSTS.') 

, 'MODEL WITH ACQUISITION COSTS; NO MOVEMENT COSTS.') 

, 'MODEL WITH MOVEMENT COSTS; NO ACQUISITION COSTS.') 

, 'THERE IS NO TOTAL ACQUISITION COST LIMIT.') 
.'INDIVIDUAL ACQUISITION COST IS $',I9,/,' A TOTAL', 

' ACQUISITION COST LIMIT OF $',I9,' IS ENFORCED.') 

, 'THERE IS NO TOTAL MOVEMENT COST LIMIT.',/) 
.'INDIVIDUAL MOVEMENT COST IS $',I9,/,' A TOTAL', 

' MOVEMENT COST LIMIT OF $',I9,' IS ENFORCED.',/) 
.'INDIVIDUAL ACQUISITION COST IS $',I9) 

.'INDIVIDUAL MOVEMENT COST IS $',I9) 



FORMAT (' 



FORMAT (' 
y 

FORMAT ( ' 



FORMAT ( ' 
h 

FORMAT( ' 
FORMAT (' 
FORMAT ( ' 
FORMAT (' 
FORMAT ( ' 
FORMAT (' 
h 

FORMAT ( ' 
FORMAT (' 



FORMAT (' 

FORMAT (' 

RETURN 

END 

SUBROUTINE OUT(NUMDIST , NUMPLAT , NUMISS , HRS , HOURS , COST , CNT , 

+ MISSIONS , DISTRICT , PLATFORM , TYPE , CATA , CATB , VALA , 

+ VALB , ATOS , STOA , QTY , PCNT , INDP , LOCINDP , WORDA , WORDB) 

********************************************************************** 
* 

* 

* 

* 

* 

* 

* 

* 

* 



* 
* 
* 
* 
* 
* 
* 
* 
* 

********************************************************************** 



SUBROUTINE 

AUTHOR 

WRITTEN 

MODIFIED 



OUTPUT DATABASE TO FILE OR SCREEN 
LT J.E. TOMKO , USCG 
12 JUNE 1991 
12 JUNE 1991 
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* PURPOSE: OUTPUT THE DATABASE IN A USER FRIENDLY FORMAT TO A FILE * 

* OR TO THE SCREEN. * 

********************************************************************** 
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INTEGER NUMDIST, NUMPLAT(4) , NUMISS, HRS(24), INDP, LOCINDP(7) 
INTEGER HOURS (7 ,10,6) , QTY(24,7), FILE 
INTEGER CNT, CATB(24,3), CATA(24,10) 

REAL COST (24) , PCNT 

REAL VALA(4) , VALB(4) , ATOS(IO), STOA(IO) 

CHARACTER MISSIONS(10)*5 , DISTRICT(7)*10 , PLATFORM (24) *8 
CHARACTER TYPE(24)*3, ANS*1, WORDA(4)*18, WORDB(4)*18 
CHARACTER CODE(3)*3, STR(3)*16 
LOGICAL ERROR 

DATA STR/'DISTRICT SURFACE' DISTRICT AIR ' , 'AREA SURFACE '/ 
DATA CODE/ 'IDS', ' IDA ' , ' IAS ' / 

ERROR - .TRUE. 

IF (CNT .GE. 1) THEN 
WRITE(* , 55) 

ENDIF 

DO 10 WHILE (ERROR) 

PRINT* ,' OUTPUT TO (F)ILE OR (S)CREEN: ' 

READ(* ,11) ANS 
IF (ANS .EQ. 'F') THEN 
IF (CNT .GE. 1) THEN 
WRITE(* ,66) 

ENDIF 

ERROR - .FALSE. 

FILE - 31 

OPEN(31 , FILE— '/DATA LISTING Al') 

CNT - CNT + 1 

ELS El F (ANS .EQ. 'S') THEN 
ERROR - .FALSE. 

FILE - 6 
ELSE 

WRITE(* , 22) ANS 
ENDIF 
ENDDO 

ERROR - .TRUE. 

DO 20 WHILE (ERROR) 

WRITE(* ,33) 

READ(*,11) ANS 
IF (ANS .EQ. 'M') THEN 

CALL OUTMISS (FILE , HOURS , MISSIONS , DISTRICT , NUMDIST , NUMISS ) 
ELSEIF (ANS .EQ. 'P') THEN 

CALL OUTPLAT (FILE , NUMPLAT , PLATFORM , TYPE , COST , QTY , CATA , 

CATB , VALA , VALB , WORDA , WORDB , STR , CODE , NUMI S S , 

MI S S IONS , NUMDI ST , DI STRI CT , HRS ) 

ELSEIF (ANS .EQ. 'D') THEN 

CALL OUTDI ST ( FILE , HOURS , MISSIONS , DI STRICT , NUMDIST , 

NUMISS , PLATFORM , QTY , TYPE , STR , CODE , INDP , 

LOCINDP, NUMPLAT) 

ELSEIF (ANS .EQ. 'T') THEN 

CALL OUTPAR ( FILE , ATOS , STOA , NUMISS , MISSIONS , VALA , VALB , WORDA , 
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20 



11 

22 

33 



44 

55 



66 



+ WORDB , PCNT ) 

ELS El F (ANS .EQ. 'A') THEN 

CALL OUTMISS (FILE, HOURS, MISSIONS, DISTRICT, NUMDIST, NUMISS) 
CALL OUTPLAT ( FILE , NUMPLAT , PLATFORM , TYPE , COST , QTY , CATA , 

+ CATB , VALA , VALB , WORD A , WORDB , S TR , CODE , NUMI S S , 

+ MI SS IONS, NUMDI ST, DISTRICT, HRS) 

CALL OUTDIST(FILE, HOURS, MISSIONS, DISTRICT, NUMDIST, 

+ NUMIS S , PLATFORM , QTY , TYPE , STR , CODE , INDP , 

+ LOCINDP, NUMPLAT) 

CALL OUTPAR ( FILE , ATOS , STOA , NUMI S S , MI S S IONS , VALA , VALB , WORDA , 

+ WORDB, PCNT) 

ELSEIF (ANS .NE. 'E') THEN 
WRITE(*,22) ANS 
ELSE 

ERROR - .FALSE. 

END IF 
ENDDO 

IF (FILE .EQ. 31) THEN 
WRITE(* ,44) 

CLOSE(31) 

END IF 
FORMAT (Al) 

FORMAT ( ' ','*** ERROR ***',/,' (',Al,') IS NOT A VALID ENTRY',/) 
FORMAT ( ' ' ,17X, 'OUTPUT MENU' ,/, 

+ ' THE FOLLOWING DATA ARE AVAILABLE FOR OUTPUT:',/, 

+ ' (D)ISTRICT DATA' ,/, ' (P)LATFORM DATA' ,/, ' (M)ISSION ’, 

+ 'DATA',/,' (T) ABLE DATA' ,/, ' (A)LL DATA' ,/, ' (E)XIT OUT' , 

+ 'PUT ROUTINE',/,' ENTER YOUR CHOICE (D,P,M,T,A, OR E) : ' ) 

FORMAT ( ' ', 'OUTPUT SENT TO FILE: DATA LISTING Al' ,/) 

FORMAT ( ' ','*** WARNING ***',/,' YOU HAVE ALREADY SENT OUTPUT ', 
+ 'TO THE FILE.',/,' SELECTING (F)ILE OUTPUT WILL ERASE ', 

+ 'THE LAST FILES PRINTED.',/) 

FORMAT (' ' , ' (F)ILE OUTPUT SELECTED',/,' OVERWRITING PREVIOUS ', 
+ 'OUTPUT.',/) 

RETURN 

END 



SUBROUTINE OUTMISS ( FILE , HOURS , MIS S IONS , DI STRI CT , NUMDIST , NUMISS ) 
********************************************************************** 



* * 

* SUBROUTINE: WRITE MISSION DATA TO SCREEN OR FILE * 

* * 

* AUTHOR : LT J.E. TOMKO, USCG * 

* * 

* WRITTEN : 12 JUNE 1991 * 

* * 

* MODIFIED : 12 JUNE 1991 * 

* * 
********************************************************************** 

* PURPOSE: WRITE THE MISSION DATA TO OUTFILE (SCREEN OR FILE) IN A * 

* READABLE FORM FOR USER. * 
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********************************************************************** 
INTEGER I, J, K, FILE, NUMDIST, NUMISS, H0URS(7 , 10 , 6) 

INTEGER TOTHRS, GRAND (7), CNT 
CHARACTER MISSIONS ( 10)*5 , DISTRICT(7)*10 
CNT - 100 

DO 10 I - 1, NUMISS 
DO 20 K - 1,7 
GRAND(K) - 0 
20 CONTINUE 

IF (CNT .GE. 63-NUMDIST) THEN 
WRITE (FILE, 11) MISSIONS(I) 

CNT - 5 
ELSE 

WRITE (FILE, 44) MISSIONS(I) 

CNT - CNT + 5 
END IF 

DO 30 J - 1, NUMDIST 
TOTHRS - 0 
DO 40 K - 1,6 

TOTHRS - TOTHRS + HOURS (J, I, K) 

40 CONTINUE 

WRITE (FILE, 22) DISTRICT(J) , (HOURS (J , I ,K) , K-1,6), TOTHRS 
CNT - CNT + 1 

GRAND ( 7 ) - GRAND(7) + TOTHRS 
DO 50 K - 1,6 

GRAND(K) - GRAND(K) + HOURS (J, I, K) 

50 CONTINUE 

30 CONTINUE 

WRITE (FILE ,33) (GRAND (K) , K-1,7) 

CNT - CNT + 2 

10 CONTINUE 

11 FORMAT('l' , 'MISSION: ’ ,A5 ,/, 12X, ' SURFACE HOURS REQUIRED' , 6X , 

+ 'AIR HOURS REQUIRED’ ,/,15X, ’CLASS’ ,21X, 'CLASS’ ,/,HX, 

+ 'DISTRICT’ , 7X , ’ 1 ’ , 6x , ' 2 ' , 6x , ' 3 ' , lOx , 

+ '1' ,7X, '2’ ,5X, '3' ,8X, 'TOTAL' ,/,lX, 

+ 8(’-’),2X,3(2X,5(’-’)) , 4X,3(2X,5(’-’)),6X,6(’-’)) 

22 FORMAT ( ’ ’ , AlO , 2X, 3 (15 , 2X) ,4X, 3(15 , 2X) , 3X, 17) 

33 FORMAT ( ' ', 12X, 5( ’-'), 2(2X, 5 ('-’)), 6X, 3(5( '-'), 2X) ,4X, 

+ 6( ’-'),/, IX, ’TOTAL’ , 6X, 3(16 , IX) , 4X, 3(16 , IX) , 3X, 18) 

44 FORMAT (/, ’ MISSION: ’ ,A5 ,/, 12X ,' SURFACE HOURS REQUIRED’ , 6X, 

+ ’AIR HOURS REQUIRED’ ,/,15X, ’CLASS' , 

+ 10X, ’CLASS’ ,/, IX, ’DISTRICT’ ,7X, 

+ '1’ ,6X, '2' ,6X, '3' ,10X, '1' ,7X, '2’ ,5X, '3’ ,8X, 'TOTAL' ,/, 

+ 1X,8('-'),2X,3(2X,5('-')),4X,3(2X,5(’-’)),6X,6(’-')) 

RETURN 
END 



SUBROUTINE OUTPLAT( FILE , NUMPLAT , PLATFORM , TYPE , COST , QTY , CATA , 

+ CATB , VALA , VALB , WORDA , WORDB , STR , CODE , NUMISS , 

+ MISSIONS, NUMDIST, DISTRICT, HRS) 

********************************************************************** 
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* 


SUBROUTINE: 


OUTPUT PLATFORM DATA 


* 


* 






* 


* 


AUTHOR : 


LT J.E. TOMKO , USCG 


* 


* 






* 


* 


WRITTEN : 


12 JUNE 1991 


★ 


* 






* 


* 


MODIFIED : 


12 JUNE 1991 


★ 


* 






* 






* PURPOSE: WRITE THE PLATFORM DATA TO OUTFILE (SCREEN OR FILE) IN A * 

* READABLE FORMAT. * 

********************************************************************** 



20 



30 



40 



50 



INTEGER I, J, K, NUMDIST, NUMPLAT(4) , NUMISS , HRS (24) 
INTEGER QTY(24, 7) , CATB(24,3), CATA(24,10), FILE, CNT 
REAL VALA(4) , VALB(4) , COST(24) 

CHARACTER MISSIONS(10)*5 , PLATFORM ( 24) *8 , STR(3)*16 
CHARACTER CODE(3)*3, TYPE(24)*3, WORDA(4)*18, WORDB(4)*18 
CHARACTER DISTRICT(7)*10 
CNT - 100 

DO 10 I - 1, NUMPLAT ( 1 ) 

DO 20 J - 1,3 

IF (TYPE(I) .EQ. CODE(J) ) K - J 
CONTINUE 

IF (CNT .GE. 55- (NUMMISS+INT(NUMDIST/4)+l) ) THEN 
WRITE(FILE.ll) PLATFORM(I) , STR(K) , COST(I) ,HRS(I) 

CNT - 5 
ELSE 

WRITE (FILE ,22) PLATFORM ( I ) ,STR(K) ,COST(I) ,HRS(I) 

CNT - CNT + 5 
ENDIF 

WRITE (FILE ,33) 

CNT - CNT + 4 
DO 30 J - 1, NUMISS 

WRITE (FILE, 44) MISSIONS(J) ,WORDA(CATA(I , J) ) , 

+ INT( (VALA(CATA(I , J) )*100 . 0)+0 . 5) 



CNT = CNT + 1 
CONTINUE 
WRITE (FILE, 55) 

CNT - CNT + 4 
DO 40 J - 1, 3 

WRITE(FILE,66) J ,WORDB(CATB(I , J) ) , 

+ INT( (VALB(CATB(I , J) )*100 . 0)+0 . 5) 



CNT - CNT + 1 
CONTINUE 

WRITE (FILE ,77) PLATFORM ( I ) 

CNT - CNT + 2 

DO 50 K - 1, (NUMDIST/3) + 1 

WRITE ( FILE ,88) (DISTRICT(J) ,QTY(I , J) , J-(K-1)*3+1, 

+ MIN((K-l)*3+3, NUMDIST)) 



CNT - CNT + 1 
CONTINUE 



10 CONTINUE 

11 FORMAT ( ' 1 ' , / , ' PLATFORM NAME: ' ,A8, 3X, 'TYPE: ' ,A16 ,/, 10X, 'YEARLY' 

+ ,' OPERATING DATA: ' ,/,6X, 'COST: ' , F6 . 2 , 7X, 'HOURS : ',14,/) 

22 FORMAT (//, ' PLATFORM NAME: ' ,A8,3X, 'TYPE: ' ,A16 ,/, 10X, 'YEARLY ', 

+ ' OPERATING DATA: ' ,/,6X, 'COST: ' .F6.2.7X, 'HOURS: ',14,/) 

33 FORMAT ( ' ', 'MISSION PROFILE DATA' ,/, IX, 20( '-'),/,' MISSION' , 

+ 4X, 'DESCRIPTION' ,9X, 'PERCENTAGE' ,/, IX, 7( ' - ' ) , 2X, 16( ' - ' ) , 

+ 6X, 10( ' - ' ) ) 

44 FORMAT ( ' ' , A5 ,4X, A18 , 6X, 13 , ' X ' ) 

55 FORMAT(/, ' CLASS PROFILE DATA' , ,/, IX, 18 ( '-'),/, ' CLASS', 7x, 

+ 'DESCRIPTION' ,9X, 'PERCENTAGE' ,/, IX, 9( ' - ' ) , 2X, 16( ' - ' ) , 

+ 6X,10('-')) 

66 FORMAT ( ' ' ,4X, II , 6X, A18 , 6X, 13 , ' X' ) 

77 FORMAT(/, ' CURRENT LOCATION OF ALL ',A8,' PLATFORMS:') 

88 FORMAT ( ' ' , 3X, 3 (A10 , 2X, 13 , 5X) ) 

RETURN 

END 



SUBROUTINE OUTD I ST ( FI LE , HOURS , MI S S I ONS , DI STRI CT , NUMDI ST , 

+ NUMISS , PLATFORM , QTY , TYPE , STR , CODE , INDP , 

+ LOCINDP .NUMPLAT) 

********************************************************************** 



* * 

* SUBROUTINE: WRITE DISTRICT DATA TO SCREEN OR FILE * 

★ * 

★ AUTHOR : LT J.E. TOMKO , USCG * 

* * 

* WRITTEN : 12 JUNE 1991 * 

★ ★ 

* MODIFIED : 12 JUNE 1991 * 

★ * 

* PURPOSE: WRITE DISTRICT DATA TO OUTFILE (SCREEN OR FILE) IN A * 

* READABLE FORMAT. * 



INTEGER I, J, K, NUMDIST, NUMPLAT (4) , NUMISS 

INTEGER HOURS (7 ,10,6), QTY(24,7), FILE, INDP, LOCINDP(7) 

INTEGER TOTHRS, GRAND (7), CNT 

CHARACTER MISSIONS ( 10)*5 , DISTRICT(7)*10, PLATFORM(24)*8 
CHARACTER CODE(3)*3, STR(3)*16, TYPE(24)*3 
CNT - 100 

DO 10 I - 1, NUMDIST 
DO 20 K - 1,7 
GRAND (K) - 0 
20 CONTINUE 

IF (CNT .GE. 59- (NUMMISS+INT(NUMPLAT(l)/4)+l) ) THEN 
WRITE (FILE, 11) DISTRICT(I) 

CNT - 5 
ELSE 

WRITE (FILE, 44) DISTRICT(I) 

CNT - CNT + 5 
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ENDIF 

DO 30 J - 1, NUMISS 
TOTHRS - 0 
DO 40 K - 1,6 

TOTHRS - TOTHRS + HOURS(I,J,K) 

40 CONTINUE 

WRITE(FILE,22) MISSIONS(J), (HOURS (I , J ,K) , K-1,6), TOTHRS 
CNT - CNT + 1 

GRAND ( 7 ) - GRAND ( 7 ) + TOTHRS 
DO 50 K - 1,6 

GRAND(K) - GRAND(K) + HOURS(I,J,K) 

50 CONTINUE 

30 CONTINUE 

WRITE (FILE ,33) (GRAND (K) , K-1,7) 

CNT - CNT + 2 

WRITE (FILE, 55) DISTRICT(I) 

CNT - CNT + 4 

DO 60 K - 1, (NUMPLAT(l)/4) + 1 

WRITE(FILE , 66) (PLATFORM(J) ,QTY(J , I) , J=(K-1)*4+1, 

+ MIN( (K-l)*4+4 ,NUMPLAT(1) ) ) 

CNT - CNT + 1 
60 CONTINUE 

10 CONTINUE 

11 FORMAT( '1' , 'DISTRICT: ' ,A10 ,/, 12X, ' SURFACE HOURS REQUIRED' , 6X, 

+ 'AIR HOURS REQUIRED' ,/,15X, 'CLASS' , 

+ 10X, 'CLASS' ,/, IX, 'MISSION' ,8X, 

+ '1' ,6X, '2' ,6X, '3' ,10X, '1' , 7X, ' 2 ' , 5X, ' 3 ' , 8X, 'TOTAL' ,/, 

+ 1X,8('-'),2X,3(2X,5('-')), 4X ,3(2X,5('-')),6X,6(’-')) 

22 FORMAT ( ' ' ,A5 , 7X, 3(15 , 2X) ,4X, 3(15 , 2X) , 3X, 17) 

33 FORMAT ( ' ', 12X, 5 ('-'), 2 (2X, 5 ('-')), 6X, 3(5 ('-'), 2X) ,4X, 

+ 6 ('-'),/, IX , 'TOTAL' , 6X, 3(16, IX), 4X, 3(16, IX) , 3X, 18) 

44 FORMAT ( / , ' DISTRICT: ' ,A10 ,/, 12X, ' SURFACE HOURS REQUIRED' , 6X, 

+ 'AIR HOURS REQUIRED ' , / , 15X , ' CLASS ' , 

+ 10X, 'CLASS' ,/, IX, 'MISSION' ,8X, 

+ '1' ,6X, '2' ,6X, '3' ,10X, '1' ,7X, '2' ,5X, '3' ,8X, 'TOTAL' ,/, 

+ 1X,8('-'),2X,3(2X,5('-')),4X,3(2X,5('-')),6X,6('-')) 

55 FORMAT ( / , ' CURRENT PLATFORMS LOCATED WITHIN ',AlO,' DISTRICT ' 

+ 'BOUNDARIES: ' ,/ , IX, 4( 'NAME' ,8X, 'QTY' , 3X) ,/,lX, 

+ 4(9('-'),3X,3('-'), 3X) ) 

66 FORMAT ( ' ' ,4 (A8 , 5X, 12 , 3X) ) 

RETURN 

END 
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SUBROUTINE OUTPAR(FILE ,ATOS , STOA , NUMISS .MISSIONS ,VALA,VALB , 

+ WORDA , WORDB , PONT ) 

********************************************************************** 



* * 

★ SUBROUTINE: WRITE PARAMETER DATA TO SCREEN OR FILE * 

★ * 

* AUTHOR : LT J . E . TOMKO, USCG * 

* * 

★ WRITTEN : 12 JUNE 1991 * 

★ * 

* MODIFIED : 12 JUNE 1991 * 

★ ★ 
********************************************************************** 

★ PURPOSE: WRITE THE PARAMETER DATA TO OUTFILE (SCREEN OR FILE) IN * 

* A READABLE FORMAT. * 



★******************************** 

INTEGER I, J, K, FILE, NUMISS 

REAL VALA(4) , VALB(4) , ATOS(IO), STOA(IO) 

CHARACTER MISSIONS(10)*5 , WORDA(4)*18, WORDB(4)*18 
WRITE (FILE, 11) (MISSIONS(I) , 1=1, NUMISS) 

WRITE (FILE, 33) (INT(ATOS(I)*100.0+0. 5) , 1-1, NUMISS) 

WRITE(FILE, 22) (MISSIONS(I) , 1-1, NUMISS) 

WRITE (FILE, 33) (INT(STOA(I)*100.0+0.5) , 1=1, NUMISS) 

WRITE (FILE, 44) 

WRITE (FILE, 66) 

DO 10 I = 1,4 

WRITE (FILE, 77) I, INT (VALA( I )*100 . 0) , WORDA(I) 

10 CONTINUE 
WRITE (FILE ,55) 

WRITE (FILE ,66) 

DO 20 I = 1,4 

WRITE (FILE ,77) I, INT(VALB(I)*100.0) , WORDB(I) 

20 CONTINUE 

11 FORMAT('l' , 'AIRCRAFT CAN CURRENTLY DO THE FOLLOWING AMOUNT',/, 

+ ' OF SURFACE RESOURCE HOURS IN EACH MISSION AREA:',/, 

+ 4X, 10(A5 , 2X) ) 

22 FORMAT (//,' CUTTERS CAN CURRENTLY DO THE FOLLOWING AMOUNT',/, 

+ ' OF AIR RESOURCE HOURS IN EACH MISSION AREA:',/, 

+ 4X, 10 (A5 , 2X) ) 

33 FORMAT ( ' ', 2X, 10(13 ,' X ', 3X) ,/) 

44 FORMAT (//, ' THE FOLLOWING VALUES AND DESCRIPTIONS ARE DEFINED ', 
+ 'FOR',/,' MISSION ASSIGNMENT CODES:') 

55 FORMAT (//, ' THE FOLLOWING VALUES AND DESCRIPTIONS ARE DEFINED ', 
+ 'FOR',/,' CLASS CODES:') 

66 FORMAT ( ' ', 5X, ' CODE' , 3X, 'VALUE' , 8X, ' DESCRIPTION' ,/, 6X,4( '-') , 

+ 3X , 5 ( ' - ' ) , 4X , 18 ( ' - ' ) ) 

77 FORMAT ( ' ' , 7X, II , 5X , 13 , ' % ' , 5X, A18) 

RETURN 

END 
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SUBROUTINE EDITPLAT (NUMDIST , NUMPLAT , NUMISS , HRS , HOURS , COST , CATA , 

+ CATB .MISSIONS , DISTRICT , PLATFORM , TYPE , QTY , 

+ WORDA , WORDB , VALA , VALB ) 

*********************************************************************** 



★ 






* 


* 


SUBROUTINE : 


EDIT PLATFORM LISTING 


* 


★ 






* 


* 


AUTHOR : 


LT J.E. TOMKO , USCG 


★ 


* 






* 


★ 


WRITTEN : 


16 MAY 1991 


* 


★ 






★ 


★ 


MODIFIED : 


30 AUGUST 1991 


* 



*********************************************************************** 

* PURPOSE: MENU DRIVEN USER INTERFACE TO ALLOW FOR MODIFICATION OF * 

* THE COAST GUARD PLATFORM DATA. USER MAY ADD, DELETE, CHANGE, AND * 

* LIST THE CURRENT PLATFORMS ON FILE. THE MAXIMUM NUMBER OF PLATFORMS * 

* WAS SET AT 24 TOTAL. THIS IS TO AVOID EXTREMELY LARGE GAMS MODELS. * 

* ROUTINE IS WRITTEN WITH ERROR CHECKING CODE TO AVOID DUPLICATION * 

* OF PLATFORM NAMES OR ATTEMPTING DELETION OF NONEXISTING NAMES. * 

*********************************************************************** 

INTEGER I, J, K, NUMPLAT (4 ) , LOCATE, FINDP , HOLD, NUMDIST 
INTEGER NUMISS, HRS (24) , HOURS(7 , 10 , 6) , QTY(24,7) 

INTEGER CATB (24,3) , CATA(24,10), PICK, COUNT, LOOP 
REAL COST(24) , VALA(4) , VALB(4) 

CHARACTER MISSIONS (10)*5 , DISTRICT(7)*10 , C0DE(4)*3 
CHARACTER PLATFORM (24) *8 , TEMP*8 , ANS*1, TYPE (24) *3, CG*26 
CHARACTER WORDA(4)*18, WORDB(4)*18 
DATA CODE/' ' , ' IDS ' , ' IDA' , ' IAS ' / 

PICK - 0 

DO 30 WHILE (PICK .EQ. 0) 

WRITE(*, 155) 

READ (*,66) ANS 
IF (ANS .EQ. 'A') THEN 

CG - 'DISTRICT SURFACE PLATFORMS' 

PICK - 2 

ELSEIF (ANS .EQ. 'B') THEN 
CG ■= 'DISTRICT AIR PLATFORMS 
PICK - 3 

ELSEIF (ANS .EQ. 'C') THEN 
CG - 'AREA SURFACE PLATFORMS 
PICK - 4 

ELSEIF (ANS .EQ. 'D') THEN 
CG - 'COAST GUARD PLATFORMS 
PICK - 1 
ELSE 

WRITE(* , 166) ANS 
ENDIF 

30 ENDDO 

WRITE (*,22) NUMPLAT ( PICK) , CG 
IF (PICK .EQ. 1) THEN 
LOOP = NUMPLAT(l)/6 
DO 10 I - 1, LOOP 
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WRITE(* , 199) (PLATFORM(J) , J-1+(I -1)*6 , 6+(I -1)*6) 

10 CONTINUE 

IF (M0D(NUMPLAT(1) ,6) .NE. 0) THEN 

WRITE(* , 199) (PLATFORM(J) , J -L00P*6+1, NUMPLAT ( 1) ) 

ENDIF 

ELSE 

COUNT - 0 
I - 1 

WRITE(*,188) 

DO 15 WHILE (COUNT .NE. NUMPLAT(PICK) ) 

IF (TYPE (I) .EQ. CODE(PICK) ) THEN 
COUNT - COUNT + 1 

WRITE(*, 177) PLATFORM(I) ,HRS(I) ,COST(I) 

ENDIF 
1 - I + 1 
15 ENDDO 

ENDIF 
ANS - 'A' 

DO 20 WHILE (ANS .NE. ’E') 

WRITE(* , 55) 

READ(*, 66) ANS 
IF (ANS .EQ. 'A') THEN 

IF (NUMPLAT(l) .EQ. 24) THEN 
WRITE(* , 144) 

ELSE 

WRITE(*,77) 24 - NUMPLAT(l) 

WRITE(*,88) 

READ(* , 99) TEMP 

LOCATE - FINDP (TEMP, PLATFORM, NUMPLAT) 

IF (LOCATE .EQ. 0) THEN 

PLATFORM (NUMPLAT(1)+1) - TEMP 
NUMPLAT (1) - NUMPLAT (1) + 1 

CALL PLAT DATA ( PLATFORM , NUMPLAT , TYPE , HRS , COST , QTY , CATA , 
+ CATB , WORDA , WORDB , VALA , VALB .MISSIONS, 

+ NUMISS.NUMDIST, DISTRICT) 

ELSE 

WRITE(* , 133) TEMP 
ENDIF 
ENDIF 

ELSEIF (ANS .EQ. 'L') THEN 
WRITE(*,22) NUMPLAT ( PI CK) , CG 
WRITE(* , 188) 

IF (PICK .EQ. 1) THEN 
DO 50 I - 1, NUMPLAT ( 1 ) 

WRITE(* , 177) PLATFORM(I) ,HRS(I) ,COST(I) 

50 CONTINUE 

ELSE 

COUNT - 0 
I - 1 

DO 25 WHILE (COUNT .NE. NUMPLAT(PICK) ) 

IF (TYPE(I) .EQ. CODE(PICK) ) THEN 
COUNT = COUNT + 1 
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WRITE(* , 177 ) PLATFORM(I) ,HRS(I) ,COST(I) 

ENDIF 
I - I + 1 
25 ENDDO 

ENDIF 

ELSEIF (ANS .EQ. 'D') THEN 
IF (NUMPLAT(l) .LT. 2) THEN 
PRINT*,'*** ERROR ***' 

PRINT* , ' AT LEAST 1 PLATFORM MUST REMAIN DEFINED ' 

ELSE 

PRINT* , ' ENTER PLATFORM NAME TO BE DELETED ' 

READ(* , 99) TEMP 

LOCATE - FINDP (TEMP, PLATFORM, NUMPLAT) 

IF (LOCATE .EQ. 0) THEN 
WRITE(* , 111) TEMP 
ELSE 

DO 45 I - 2,4 

IF (TYPE(LOCATE) .EQ. CODE(I)) THEN 
NUMPLAT ( I ) - NUMPLAT ( I ) - 1 
ENDIF 

45 CONTINUE 

IF (LOCATE .NE. NUMPLAT ( 1 ) ) THEN 
DO 40 I - LOCATE, NUMPLAT (1)-1 
PLATFORM(I) - PLATFORM(I+l) 

TYPE(I) - TYPE(I+1) 

HRS(I) - HRS(I+1) 

COST ( I ) - COST (1+1) 

DO 41 J - 1,3 

CATB(I , J) - CATB(I+1 , J) 

41 CONTINUE 

DO 42 J - l.NUMISS 

CATA ( I , J ) -= CATA( 1+1 , J ) 

42 CONTINUE 

40 CONTINUE 

ENDIF 

NUMPLAT ( 1 ) - NUMPLAT ( 1 ) - 1 
ENDIF 
ENDIF 

ELSEIF (ANS .EQ. 'C') THEN 

PRINT*, 'ENTER NAME OF PLATFORM TO BE CHANGED' 

READ (*,99) TEMP 

LOCATE - FINDP (TEMP, PLATFORM, NUMPLAT) 

IF (LOCATE .EQ. 0) THEN 
WRITE(* , 111) TEMP 
ELSE 

CALL CHGPLAT (NUMDIST , NUMPLAT , NUMISS , MISSIONS , DISTRICT , 
+ PLATFORM , CATA , CATB , VALA , VALB , WORDA , WORDB , 

+ HRS, QTY, COST, LOCATE, TEMP) 

ENDIF 

ELSEIF (ANS .EQ. 'E') THEN 

PRINT*, 'EXITING EDIT PLATFORM ROUTINE' 

ELSE 
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